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THE CATALYTIC REMOVAL OF AMMONIA AND NITROGEN 

OXIDES FROM SPACECABIN ATMOSPHERES 

By A. J. Gully, R. R. Graham, J. E. Halligan, 
and P. C. Bentsen 

Chemical Engineering Department 
Texas Tech University 


1 .0 SUMMARY 


Investigations were made on methods for the removal of ammonia and to 
a lesser extent nitrogen oxides in low concentrations from air. The cata- 
lytic oxidation of ammonia was studied over a temperature range of 250°F 
to 600°F and a concentration range from 20 ppm to 500 ppm. Of the catalysts 
studied, 0.5% ruthenium supported on alumina was found to be superior. This 
material is active at temperatures as low as 250°F and was found to produce 
much less nitrous oxide than the other two active catalysts, platinum on 
alumina and Hopcalite. 

A quantitative design model was developed which will permit the per- 
formance of an oxidizer to be calculated. The ruthenium was found to be 
relatively insensitive to low concentrations of water and to oxygen con- 
centration between 21% and 100%. Hydrogen sulfide was found to be a poison 
when injected in relatively large quantities. The ruthenium catalyst was 
found to have a very stable activity under normal operating conditions with 
one charge having been used for several weeks without a loss in activity. 

The adsorption of ammonia by copper sulfate treated silica gel was in- 
vestigated at temperatures of 72°F and 100°F. Ammonia feed concentrations 
investigated ranged from 80 ppm to slightly over 500 ppm. A quantitative 
model was developed for predicting adsorption bed behavior. Limited investi- 
gation indicated that the treated silica gel could be thermally desorbed with 
little loss of capacity. 


2.0 INTRODUCTION 


In a closed system such as a spacecabin, the volume of free atmosphere 
per inhabitant is extremely small resulting in large increases in concentra- 
tion of contaminants with the addition of very small amounts of contaminant 
gases. If the contaminant is metabolically produced (as is ammonia) then 
very reliable systems for contaminant removal must be included in the con- 
taminant control system. The release of nitrogen containing compounds into 
the atmosphere of spacecabins is an unavoidable occurrence even though con- 
siderable effort is made to control these compounds at the source of the re- 
lease. Once the nitrogenous compounds are released into the atmosphere of 
the spacecabin, their concentration must be controlled by the trace contam- 
inant removal system in order to maintain a life supporting environment. 

The contaminant of particular concern in this investigation is ammonia. Am- 
monia production rates are estimated to be as high as 1.2 x 1CT 4 Ib/hr. with- 
in a 2000 cubic foot cabin volume. 

Two different control methods were evaluated in this investigation. 

These methods were the catalytic oxidation of ammonia and the adsorption of 
ammonia on treated silica gel. The removal of ammonia by catalytic oxidation 
to yield nitrogen and water is theoretically a very attractive method of 
controlling this contaminant provided no undersirable reaction products are 
produced. Currently proposed trace contaminant control systems contain a 
catalytic burner through which the cabin air is circulated. The products of 
ammonia oxidation within the catalytic burner could include nitrogen dioxide 
(N0 2 ), nitrous oxide (N 2 0), and/or nitric oxide (NO). The particular nitro- 
gen oxides produced and their rates of production undoubtedly depend on the 
type of catalyst used and on the temperature at which the oxidation is carried 
out. For this reason the reaction rates and product distribution for the cat- 
alytic oxidation of ammonia were explored over a rather wide range of operating 
conditions and catalytic materials. 

The second method of ammonia removal investigated, adsorption on copper 
sulfate treated silica gel, is also attractive provided the capacity of the 
solid for ammonia is sufficiently high and provided regeneration of the sat- 
urated bed can be made rapidly and without significant loss of capacity. 

The oxides of nitrogen and ammonia are all undesirable contaminants in 
a life supporting atmosphere. The presence of nitrogen dioxide in particular 
must be minimized. The maximum allowable concentration levels measured at 
one atmosphere and 72°F for continuous exposure to these materials are reported 
(ref. 1) as follows: 


NH 3 5.0 ppm by volume 

N 2 0 26.0 ppm by volume 

NO 26.0 ppm by volume 

N0 2 0.5 ppm by volume 
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2.1 Objectives 


The objective of the work reported here was to investigate the removal 
of ammonia and its oxidation products from spacecabin or other closed system 
atmospheres. In order to pursue this objective, the following research and 
development program was followed: 

1. Analytical techniques were perfected which would permit the con- 
venient and accurate quantitative analysis of mixtures of ammonia, 
the nitrogen oxides and air. 

2. Several noble metal and metal oxide catalysts were screened in order 
to determine their activity for atmonia decomposition. 

3. The three most promising catalytic materials (ruthenium supported 
on alumina, platinum supported on alumina, and Hopcalite) were 
investigated for ammonia decomposition activity and product distri- 
bution over a wide range of conditions. 

4. The ammonia oxidation data from Step 3 was correlated in order to 
permit prediction of reactor performance under conditions expected 
within the contaminant control system of a spacecabin. 

5. The adsorption of ammonia on copper sulfate impregnated silica gel 
was investigated. A model was developed which would permit the 
prediction of adsorber performance. 


2.2 Research Outline 

The search for an effective ammonia removal technique was pursued along 
two fundamentally different routes; these were (a) catalytic oxidation at 
elevated temperature, and (b) adsorption at ambient temperature. These two 
efforts were carried out simultaneously and essentially independently. 

From very early in the investigation it was known that the major problem 
which would have to be overcome in order to make catalytic ammonia oxidation 
a viable method for ammonia removal was to determine ways of minimizing the 
production of nitrogen oxides. The two major manipulatable variables available 
for controlling the nitrogen oxides production were catalyst type and oxida- 
tion temperature. The experimental program then was to determine the general 
characteristics of several potentially active catalysts and from these try 
to pick two or three for more detailed study. Historically the most active 
ammonia oxidation catalysts known were supported noble metal. 

The optimization of reactor temperature is a crucial item. Higher 
reaction temperatures tend to produce proportionally more nitrogen oxides, 
particularly nitric oxide and nitrogen dioxide. Since nitric oxide is known 
to be easily oxidized to nitrogen dioxide at lower temperatures, these two 
oxides can be considered together. The other undersirable product, nitrous 
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oxide, is produced in most abundance at temperatures below those which yield 
nitrogen dioxide. Thus the essential problem was to determine the rates of 
ammonia removal and nitrogen oxides production over the temperature range 
potentially available for use in a catalytic removal device. 

In the ammonia adsorption phase of the project, the experimental pro- 
gram was designed to provide a series of ammonia breakthrough curves at 
constant temperature and humidity level. From these curves an isotherm could 
be determined. Breakthrough curves would also supply the necessary information 
for development of quantitative relationships for the design of an ammonia 
adsorption system. 
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3.0 SYMBOLS 



'Ai 

'Ao 

'Bo 

: N0 2 

"NoO 


D 

°K 

d p 

E' 

E f 

E NH 3 

e n 2 o 

F i 

g m 

h 

K 


k f 

ki 

(ki) 


obs 


experimental constant in Langmuir equation (ft 3 /lb) 

external area of adsorbent, ft 2 /! b 

pre-exponential term in ammonia rate expression 

pre-exponential term in N 2 0 rate expression 

external area of catalyst per unit mass (cm 2 /g cat) 
experimental constant in Langmuir equation (ft 3 /mol) 

concentration of adsorbate in bulk gas phase (mols/ft 3 ) 
dimensionless constant in BET equation 

concentration of adsorbate in pore (mols/ft 3 ) 

adsorbate concentration in equilibrium with solid surface 
(mols/ft 3 ) 

ammonia concentration in bulk phase (ppm) 
inlet concentration of A (mols A/total mols) 

outlet concentration of A (mols A/ total mols) 

outlet concentration of B (mols B/total mols) 

outlet N0 2 concentration (ppm) 

outlet N 2 0 concentration (ppm) 

ammonia concentration at gas-solid interface (ppm) 

concentration of component j 

diffusivity (cm 2 /sec) 

Knudsen diffusivity (cm 2 /sec) 

catalyst particle diameter (cm) 

constant in Arrhenius equation, activation energy (Btu/mol) 
effectiveness factor 

activation energy for ammonia reaction (Btu/lb mol/°R) 

activation energy for N 2 0 formation (Btu/lb mol/°R) 

inlet flow rate (g mol /hr) 

superficial velocity (mol /hr/ cm 2 ) 

Thiele modulus, dimensionless 

bulk gas phase mass transfer coefficient (ft/hr) 

reaction rate constant 

mass transfer coefficient (g mol/hr/atm/cm 2 ) 
true first order rate constant 

ammonia reaction rate constant 
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%0 


Re 


M 

N 

N 

n 

P 

P T 

p g 

p i 

P o 

R 

ft 


r 

r. 


NH 3 

r N 2 0 

r N0 2 

r eff 


J X 

T 

T, T' 
V 

V„ 


W 

M 

W E 

W M 

YK 


Z 

Ci 

3 

e B 

Ep 


N 2 0 formation rate constant (g mol N 2 0/hr/g cat/(ppm NH 3 ) 0 - 62 ) 

experimental rate constant 

molecular weight of diffusing species 
rate of mass transfer (g mol/hr/g cat) 

Reynolds number, dimensionless 

empirically determined reaction order 
partial pressure 
total pressure (atm) 

partial pressure in bulk phase (atm) 

partial pressure at gas solid interface (atm) 

vapor pressure of adsorbate, in. Hg 

bed radius, ft 

adsorption rate (mols gas/ft 3 solid - hr) 
pore radius (cm) 

universal gas constant (Btu/mol °R) 
reaction rate 

N 2 formation rate (g mol/hr/g cat) 
ammonia reaction rate (g mol/hr/g cat) 

N 2 0 formation rate (g mol/hr/g cat) 

N0 2 formation rate (g mol/hr/g cat) 

reaction rate in a porous catalyst when there is no internal 

diffusion resistance (g mol/hr/g cat) 

external area of a single catalyst particle (cm 2 ) 

time (hr) 

absolute temperature (°R) 

gas phase superficial velocity (ft/hr) 

total volume of a single catalyst particle (cm 3 ) 

pore volume per unit mass of catalyst (cm 3 /g cat) 

catalyst weight (g) 

adsorbate loading on solid (mols adsorbate/lb adsorbent) 
equilibrium surface concentration (mols/lb) 

monomolecular coverage of adsorbent by gas (mols/lb) 

adsorption rate constant for surface adsorption rate controlling 
the adsorption (lb/ft 3 hr) 
bed length (ft) 

empirically determined ammonia reaction order, dimensionless 
empirically determined N 2 0 formation order, dimensionless 
adsorbent void fraction, diemns ionless 

internal void fraction, dimensionless 
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r = pore diffusivity, ft 2 /hr 

y = viscosity (g/cm/sec) 

p = density (g/cm 3 ) 

p = apparent adsorbent density, lbs/ft 3 
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4.0 PERFORMANCE OF AN AMMONIA OXIDIZER 


A major objective of the present study was to develop the data and 
techniques necessary to design and predict the operating performance of a 
catalytic oxidizer for ammonia. It was realized early in the project that 
the variables of primary concern were catalyst type and operating tempera- 
ture. The activity of catalytic materials for ammonia oxidation varies 
greatly from material to material. The activity and product distribution 
of a catalyst are strong functions of reaction temperature. In the follow- 
ing sections the experimental behavior of ammonia over several catalysts is 
presented. 


4.1 Literature Survey 

A considerable amount of information regarding the high temperature 
oxidation of ammonia (above 450°C) is available due to the commercial signi- 
ficance of this reaction in the manufacture of nitric acid (ref. 2). At a 
temperature of 800°C ammonia can be oxidized almost exclusively to NO using 
a platinum catalyst. 

Much less information is available on the oxidation of ammonia at 
temperatures below 450°C. Griffiths et al_. (ref. 3) performed an infrared 
study of the adsorption and oxidation of ammonia on silica-supported platinum 
at temperatures up to 300°C. The authors operated with a 12:1 molar ratio of 
oxygen to ammonia (total pressure 10 torr) which they state is the optimum 
condition for a maximum N 2 0 yield. Up to 178°C they observed the same spec- 
tra for ammonia and oxygen as they did for ammonia alone. At 178°C and 191°C 
they observed a doublet corresponding to gaseous N 2 0. The appearance of NO 
or N0 2 was not reported. The appearance of N 2 was not reported and it is 
not known whether N 2 would have been detectable with their apparatus. 

Amano and Taylor (ref. 4) studied the decomposition of ammonia admixed 
with N 2 and H 2 only over ruthenium, palladium, and rhodium catalysts supported 
on alumina. Their temperatures ranged from 340°C and 400°C and the maximum 
ammonia decomposition observed was 10%. Hydrogen was found to inhibit the 
reaction over palladium and ruthenium while N 2 was found to be without influ- 
ence. An observed activation energy of 30 kcal/mol was reported. 

Schriber and Parravano (ref. 5) investigated the kinetics of ammonia 
oxidation over a 0.5% ruthenium catalyst supported on aluminum oxide using a 
fixed bed recirculation reactor. Their investigation was conducted at six 
temperatures in the range of 246°C-345°C and at total pressures ranging to 4 
atm. The feed composition ranged from 5.3 to 16.5 mol percent ammonia with 
approximately 10 mol percent 0 2 , 3 mol percent H 2 0, and helium as the diluent. 
The authors found that the oxidation products were N 2 , H 2 0, and N 2 0. The 
production of N0 2 was not reported and it is questionable whether it could 
have been detected. 
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Schriber and Parravano state that the relationship between N 2 , N 2 0 
and ammonia can be described by the following reaction scheme: 



Under the experimental conditions described the authors found that k 3 was 
much less than k 2 and that the oxidation could thus be represented by two 
parallel reactions occurring on the ruthenium surface. External mass trans- 
fer was not found to be a controlling factor in the kinetics of the reaction, 
and the authors reasoned that pore diffusion could not be important since 
the ruthenium metal was deposited on the exterior of the catalyst pellet. 

The following empirical rate equations were derived from the experimental 
data: 


, , 0.46 " 0 .40 

= ^ k o n 2 p nh 3 p o 2 p h 2 o 

(4.1) 

„ 1 .35 0 .46 - 0 .40 

f k 0^N 2 0 P NH 3 P 0 2 P H 2 0 

(4.2) 


where k Q = experimental rate constant 
P u = partial pressure 
r = reaction rate 

The apparent activation energies for the formation of N 2 and N 2 0 were found 
experimentally to be 12,600 and 35,000 cal /mol respectively. The authors 
were able to reconcile their experimental k Q values with theoretical values 
calculated from absolute rate theory. To u make this calculation they 
assumed that the reaction involved imide (NH) and nitroxyl (HN0) surface 
intermediates. These two surface compounds had been hypothesized in earlier 
work (ref. 2). 

The catalytic oxidation of NH 3 by air over a promoted nickel oxide 
catalyst was investigated by N. Giordano, et al^. (ref. 6). In the temperature 
range between 300° and 450°C, the reaction was reported to proceed mainly 
to N 2 and N 2 0. As more sodium or lithium promotor was added both activity 
and yield of nitrous oxide increased. Doping with tri valent ions tended to 
depress both activity and yield of nitrous oxide. 

The minimum N 2 0 yield reported at 352°C was approximately 12% and was 
obtained using NiO doped with 2% Fe 2 0 3 . However the reaction rate was also 
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relatively low as only 4% conversion of a 0.3 gmole/hr feed stream was ob- 
tained using 0.85 g of catalyst. At 441 °C the minimum N 2 0 yield was approxi- 
mately 17% using the same catalyst and flow conditions. The sodium and 
lithium doped catalysts gave 50% to 95% yields of N 2 0. 


4.2 Modelling Techniques 

Generally speaking, the two most common methods of correlating catalytic 
reaction rate data are the overall fractional order rate expressions and the 
Langmuir-Hinshelwood rate expressions. The Langmuir-Hinshelwood method at- 
tempts to describe the overall reaction in terms of the elementary processes 
that occur at or near the catalyst surface. The overall rate expressions, on 
the other hand, neglect the elementary steps and rely instead on a simple rate 
expression which contains empirical constants having little physical signifi- 
cance. An often used overall rate expression is of the form 

r = kc*j (4.3) 

J 


where r 

= reacti on ratfe 

c i 

k J ,n 

= concentration of component j 
= empirically determined constants. 


In order to reflect the influence of temperature on the reaction rate, the 
Arrhenius equation is used for the forward rate constant. 

k = Ae" E/RT (4.4) 


Kinetic data over a relatively wide range of temperature can usually be re- 
presented by the use of this relationship. 

A primary advantage of a rate expression of the form 

r = Ae _E/RT C n (4.5) 


is that the effect of the two most important variables is correlated with a 
minimum number of experimentally determined parameters. This fact is impor- 
tant since the amount of kinetic data generally available is not sufficient 
to confidently fit a complicated multi -parameter model. 

All fluid solid catalytic reactions can be broken down into the follow- 
ing series of steps which occur either in series or in parallel near the 
catalyst surface: 

1. Transport of reactants from the bulk fluid phase to the fluid-solid 
i nterface 
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2. Transport of reactants into the catalyst pores 

3. Adsorption of reactants on the catalyst surface 

4. Surface reaction 

5. Desorption of products from the catalyst surface 

6. Diffusion of products out of the catalyst pores 

7. Transport of products back into the bulk fluid phase. 

All of the above listed steps except 2 and 3 occur in series. Steps 2 and 
3 occur in parallel because adsorption and reaction occur continually as 
the reactants pass through the catalyst pores. Steps 1 and 2 are of special 
significance to this study and will be discussed in more detail. Steps 3 
through 5 will be correlated by use of the previously discussed fractional 
order rate expression. 

4.2.1 External Diffusion . - In any catalysis study it is important to 
ascertain whether or not the reaction rate is strongly influenced by the 
rate of reactant diffusion to the catalyst surface. When reactant diffusion 
to the catalyst surface is a rate limiting step then the true nature of the 
catalytic reaction, including the reaction product distribution, can be mask- 
ed if the effects of mass transport are not properly accounted for. 

Fluid passing over the surface of a catalyst develops a boundary layer 
external to the catalyst pellet in which the fluid velocity parallel to the 
surface varies from zero, at the surface, to the bulk stream velocity over 
a short distance. In the bulk stream, mixing is primarily by turbulent 
transport, whereas in the boundary layer region transport is by molecular 
diffusion. Mathematically, the rate of mass transfer from the bulk phase 
to the catalyst surface can be described by the expression 


N = k f a( p g- p i } - (4.6) 

For reactions which are occurring at steady state the rate of mass transfer 
is also equal to the rate of reaction; therefore 


r ■ k' a(P g -P 1 ) 


(4.7) 


where 


k 1 
K f 

a 

p g 

p i 


rate of mass transfer 
mass transfer coefficient 

external area of catalyst particle per unit mass 
partial pressure in bulk phase 

partial pressure at gas-solid interface. 
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The mass transfer coefficient is a function of the physical properties of 
the compound which is diffusing, as well as the nature of the flow around 
the catalyst surface. Hougen and Wilkie (ref. 7) developed the following corre- 
lations of these variables for packed beds: 


and 


J D = 


kf Pj u 2/3 -0.41 

<§y> - °- 989 <v 


M 


J D = 


'ill (i L.) 2/Z = 1 82 fN 
G m W 1,82 (N Re) 


-0 . 51 


where 


d G m 
N = P M 

N Re p 
d 


P 

M 

m 

P 

D 


= diameter of catalyst particle 
= superficial molar velocity 
= viscosity of mixture 
= mean molecular weight 
= density of reaction mixture 
= diffusivity 

= dimensionless mass transfer number. 


(4.8) 


(4.9) 


Equation (4.8) is valid for Reynolds numbers greater than 350, while Equation 
(4.9) is valid for lower values. 


Using Equations (4.8) and (4.9) it is possible to estimate the value 
of k|. Once k^. is determined and experimental rate data are available it is 

possible to estimate the concentration gradient external to the catalyst 
surface through the use of Equation (4.7). If the value of P^-P^ in Equation 

(4.7) is a small fraction of P g then the resistance of the boundary layer can 

be neglected. However if not, then the reactant partial pressure at the gas- 
solid interface must be used in Equation (4.5). 

4.2.2 Internal Diffusion . - For a given concentration driving force the 
rate of surface adsorption in heterogeneous catalysis may be high in compari- 
son with the rate of diffusion, resulting in concentration gradients down the 
catalyst pores. For this reason the interior of a catalyst pellet may be 
less effective than the exterior. 
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Diffusion in catalyst pores may occur by ordinary diffusion, Knudsen 
diffusion, or surface diffusion. Ordinary diffusion occurs if the size 
of the pore is large with respect to the mean free path of the molecules. 

If the size of the pore is small in comparison to the mean free path, then 
collisions with the pore walls become important. Most catalysts of indus- 
trial significance have pore sizes such that Knudsen diffusion predominates 
(ref. 8). The Knudsen diffusivity can be calculated from the expression 


where 



D k = Knudsen diffusivity 
Rp = pore radius 
T = absolute temperature 


M = molecular weight of diffusing species. 


(4.10) 


Surface diffusion occurs as a result of the limited mobility which adsorbed 
molecules have. Considerably less is known about surface diffusion than 
Knudsen and ordinary diffusion. Satterfield (ref. 9) states that surface dif- 
fusion cannot be important unless very weak surface adsorption bonds are 
formed . 


The effect of pore diffusion on the apparent activation energy can be 
important, especially at high temperatures. This is because the surface pro- 
cesses (adsorption and surface reaction) increase in rate exponentially with 
temperature while diffusion rates are not affected as greatly. Thus rate 
measurements at high temperatures reflect diffusional resistances much more 
than do low temperature measurements. Satterfield (ref. 9) states that the 
effect of temperature on diffusion rates is equivalent to an activation energy 
of 1 to 3 kcal/mol , and that the apparent activation energy when pore dif- 
fusion is important will be approximately one-half that of the true activa- 

Designating r ^ as the reaction rate in a porous catalyst when there is 

no internal diffusion resistance, an effectiveness factor, E f , can be defined 
as T 



(4.11) 


It can be seen from Equation (4.11) that the effectiveness factor is the ratio 
of the observed reaction rate to the reaction rate which would be observed if 
all the active catalyst surface were on the periphery of the catalyst particle. 
The effectiveness factor can also be related to the Thiele modulus (ref. 10) 


(4.12) 


r _ tanh (h) 
L f h 
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where the Thiele modulus, h, is defined as 


h " * T - / TV — IT 
2 / g k 


(4.13) 


where V p = total volume of a single catalyst particle 
= external area of a single catalyst particle 
k x = true first order reaction rate constant 
V g = pore volume per unit mass of catalyst. 

Using the above equations along with experimentally determined rates, a 
simple trial and error procedure can be used to obtain E f . 

1. A value of is assumed, and h is calculated from Equation (4.12) 

2. Knowing V p , S x> and D^, k x is calculated from Equation (4.13) 

3. A new is calculated from 


r _ ^ kl W 

where (k : ) Qbs - observed first order rate constant. 

4. If the E f calculated agrees with the E f assumed, then the correct 
value of E f has been found. If it does not. Step 1 is repeated. 

4.2.3 Limitations on Fractional Order Models . - Fractional order rate 
expressions (as with all empirical expressions) must be used with extreme 
caution outside of the range of the independent variables used in determining 
the activation energy and concentration dependency. At higher reagent con- 
centrations the value of n in Equation (4.3 ) is generally less than one. 
Since most reactions approach a dependency of one as reagent concentration 
goes to zero, it is recommended that the value of n be made one for extra- 
polations of reaction rate at concentrations much below those used in ex- 
perimentally determining n. This will almost surely result in predicted 
reaction rates which are slower than physically occur and thus an over de- 
signed reactor. However due to the uncertainty involved this conservative 
approach is warranted. 
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4.3 Preliminary Screening of Catalyst 

During the first half of the project, ammonia oxidation data was 
obtained using thirteen different catalyst materials. Shown in Table 4.1 
are the catalysts which were included in this preliminary screen matrix. 

All of the materials tested showed some activity for the removal of am- 
monia; however, platinum, ruthenium and Hopcalite seemed most promising 
and were selected for further detailed study. Each of these materials 
will be separately discussed later. 

Both the fixed-bed reactor and the stirred-tank reactor were used 
in the preliminary catalyst screening tests. In order to make a relative- 
ly easily interpretable comparison of the low temperature ammonia oxida- 
tion activity of these catalysts, a first order reaction rate approximation 
was used and the resulting rate constant was plotted versus reciprocal tem- 
perature. 

Although this procedure may not fully compensate for concentration 
effects, the procedure does permit the ranking of catalyst in terms of 
their activity. Shown in Figure 4.1 are the results of the preliminary 
screen tests. These tests were made using air with less than 500 ppm 
water. A high first order rate constant at a low temperature is of 
course the most desirable situation. As previously mentioned ruthenium, 
Hopcalite, and platinum appear to be the most active catalyst. The two 
cobalt containing catalysts (Co 0901 and Ni 1601) were both relatively 
active, however they both tend to produce a relatively large fraction of 
highly undesirable nitrogen dioxide product. "Gatalyst" , copper oxide, and 
nickel oxide all were considerably less active than the better materials. 
The iron oxide, silver and manganese dioxide were not effective showing 
only a small amount of reaction at the highest temperatures used. 


4.4 Experimental Data - Ruthenium Catalyst 

4.4.1 Ammonia Removal Kinetics . - Using the backmix reactor system 
it was possible to collect data from which the overall ammonia reaction 
rate, N 2 0 formation rate, and N0 2 formation rate could easily be calculated. 
Table 4.2 presents the raw data for reactor inlet and outlet compositions 
at five temperatures from 240°F to 442°F. The ammonia concentration in the 
reactor outlet varied between 13 and 16 ppm while the ammonia removal rate 
varied from 16 X 10" 6 to 197 X 10" 6 gmole NH 3 /hr/g catalyst between 242°F 
and 442°F. The N 2 0 formation rate was found to be approximately one order 
of magnitude lower than the total ammonia removal rate while the N0 2 forma- 
tion rate was in most cases less than half the N 2 0 rate. 

Figures 4.2 and 4.3 show the total ammonia reaction rate versus bulk 
concentration at five temperature levels. The curved lines through the 
data on the figures were obtained by a linear least squares fit to the 
three parameter equation 


In r 


nh 3 


In A mu + a Inc. 
NH 3 i 


E NH 3 

RT 


(4.14) 
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TABLE 4.1 


CATALYSTS INCLUDED IN SCREENING TESTS 


CATALYST MATERIAL 

COMPANY DESIGNATION 

SURFACE AREA m 2 /g 

2% Mn0 2 supported on silica 

Harshaw Mn-0501 

258 

0.3% Pd supported on high 
activity alumina 

Harshaw Pd-0501 

186 

0.5% Pt supported on alumina 

Englehard 0.5% Pt 
supported on alumina 

93 

0.5% Ru supported on alumina 

Englehard 0.5% Pt 
supported on alumina 

74 

4% nickel oxide, 4% cobalt 
oxide, 4% iron oxide 
supported on activated alumina 

Harshaw Ni-1601 

78 

14% nickel oxide supported on 
high activity alumina 

Harshaw Ni-0707 

140 

10% CuO supported on high 
activity alumina 

Harshaw Cu-0803 

137 

19% Mn0 2 supported on 
activated alumina 

Harshaw Mn-0201 

69 

20% Fe 2 0 3 supported on 
activated alumina 

Harshaw Fe-0301 

41 

3.5% Ag supported on inert 
al umi na 

Harshaw Ag-0101 

1 

5% Cobalt oxide, 5% copper 
oxide supported on activated 
alumina 

Harshaw Co-0901 

59 

Not reported 

General American 

Transportation 

"Gatalyst” 

Not reported 

Mixed oxides of Manganese 
and copper 

Mine Safety Appliance 
"Hopcalite" 

156 
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** not measured 
*** humidified run 
**** pure 0 2 run 



Ruthenium Catalyst 
Least Squares Fit 



Figure 4.2 - Differential Data Showing Ammonia 
Removal Rate Versus Concentration 











Where r NH 3 
a nh 3 
e nh 3 

R 

T 

c i 

a 


ammonia reaction rate 

pre-exponential term in Arrhenius expression 
activation energy for anmonia reaction 
universal gas constant 
absolute temperature 

ammonia concentration at gas-solid interface 
reaction order. 


The equation form used results in curves passing through the origin (as 
would be expected theoretically); however, caution should be used in extra- 
polation to concentrations outside of the data set. (See also Section 4.2.2). 
Although the concentration on the absicissas of Figure 4.2 and 4.3 refer to 
bulk ammonia concentration, interfacial values were used in determining the 
constants A^ , , and a. The best values of A^ H , E^ , and a were found 

by the least squares procedure to be 


(gmole NH 3 ) 

a - ? ?fl4 x in -3 

(hr ) (g cat)(ppm NH 3 )°* 69 


E 


NH 3 


10,100 


Btu 
lb mol 


a = 0.69 


The average deviation between the predicted and experimental points of 
Figures 4.2 and 4.3 was found to be 14.5%. Inherent in using an equation 
such as (4.14) is the assumption that the order of the reaction does not 
change over the temperature range of interest. This is believed to be a 
reasonable assumption in view of the relatively low temperatures and short 
temperature range involved in this investigation. 

The rate expression and experimentally determined kinetic parameters 
discussed above were used to predict the ammonia removal for an integral 
reactor. The fixed bed reactor was then operated to provide experimental 
integral reactor data to compare with the predicted ammonia removal. In 
order to make the integral reactor calculations a numerical integration 
program was utilized which included the effects of external diffusion as 
discussed in Section 4.2.1. 
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The comparison of experimental and calculated ammonia removal is shown 
in Figures 4.4 and 4.5 for reaction temperatures of 300°F and 400°F. The 
agreement shown in these tests indicates that the performance of an integral 
reactor using ruthenium on alumina catalyst can be predicted with confidence 
from the kinetic constants given above. 


4.4.2 Nitrous Oxide Production Kinetics . - Figures 4.6 through 4.9 
show the rates of formation of N 2 0 as a function of the bulk anmonia concen- 
tration at five temperatures. The curved lines shown were obtained by a 
linear least squares fit to the equation 


e n 2 o 

r N 2 0 “ A N 2 0 exp “RT - ^ c i 


(4.15) 


where r^ 


N 2 0 formation rate 


A N 2 0 = P re " ex P onent1 term 10 Arrhenius expression 
E n q = activation energy for N 2 0 formation 


3 = reaction order. 


The best values of A 


N 2 0’ 


E N 2 0’ 


and 3 were found to be 


(gtnole N 2 0) 

A M n = 0 . 1 26 

(hr)(g cat)(ppm NH 3 )° 70 



20 ,640 


Btu 
lb mol 


3 = .70 . 

The average deviation between the predicted and experimental points for N 2 0 
formation was found to be 30%. The low temperature point contributed the 
largest amount to the average deviation due to the very small reaction rates 
involved. 
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Figure 4.5 - Comparison of Calculated and Experimental 

Ammonia Conversion for Fixed Bed Reactor 
at 400°F 
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Figure 4.7 - Nitrous Oxide Production Rate 
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Just as with the ammonia removal data, the nitrous oxide production for 
an integral reactor was predicted based on the above kinetic constants. The 
fixed bed reactor was used to generate experimental nitrous oxide production 
data for comparison with the prediction. Shown in Figure 4.10 are the pre- 
dicted and experimental nitrous oxide production for temperatures of 300°F 
and 400°F. The experimental nitrous oxide yield was less at both temperatures 
than the model prediction. Thus the model may be expected to predict more 
nitrous oxide than will be actually experienced. 

Since the ammonia concentration dependency (or reaction order) of both 
the ammonia removal and nitrous oxide production are approximately the same, 
the percentage conversion to N 2 0 would not be dependent on ammonia concen- 
tration. However since the activation energy for nitrous oxide production 
is approximately twice that for ammonia removal, the percentage yield of 
nitrous oxide rises rapidly with increased temperature. 

4.4.3 Nitrogen Dioxide Production Kinetics . - Nitrogen dioxide pro- 
duction rates were found to be a strong function of reaction temperature with 
an activation energy of approximately 28,000 Btu/mole. At all temperatures 
below 400°F the conversion to nitrogen dioxide was very small. At 400°F less 
than 2% of the reacted ammonia was converted to nitrogen dioxide and at tem- 
peratures of 300°F and below less than 1% conversion of ammonia to nitrogen 
dioxide was found. 

4.4.4 Effect of Mater . - The triangular data points on the 400°F data 

of Figures 4.3 and 4.8 represent runs in which the inlet air stream to the 
reactor was humidified to 0.7 mol percent H 2 0 (23% relative humidity at 

72°F) and fed to the reactor. As can be seen from the positions of the data 

points in Figure 4.3, the overall NH 3 removal rate is not significantly de- 

pressed by the presence of H 2 0 at this concentration. Referral to Figure. 

4.8 shows that H 2 0 also exhibits no inhibitory effect on N 2 0 formation. 

These facts are not in agreement with the findings of Schriber (ref. 5) who found 
H 2 0 to be a strong inhibitor as evidenced by his proposed rate equations 


^ k o^N 2 P NH 3 


p 0* 46 p “0.40 

k o 2 k h 2 o 


and 


r N?0 “ ^NoO P NH, 1 ' 35 


p 0 .46 p “ 0. 4 0 

k o 2 k n 2 o 


4.4.5 Effect of Oxygen . - The square data points of Figures 4.3 and 4.8 
represent data taken with pure 0 2 substituted for air in the reactor inlet. 

As can be seen from Figure 4.3, the additional 0 2 does not seem to have an 
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effect on the overall ammonia removal rate. Referral to Figure 4.8 shows 
that the rate of formation of N 2 0 appears to be slightly increased due to 
the added 0 2 . If the rate of N 2 and N0 2 formation are unaffected, however, 
the additional N 2 0 formation would not be sufficient to cause a noticeable 
increase in the overall ammonia removal rate. 

4.4.6 Internal Diffusion Effects . - In order to ascertain whether or 
not pore diffusion was a significant resistance, an effectiveness factor 
was calculated using the simplified scheme presented by Smith (ref. 8) and ex- 
plained in section 4.2.2. Using an approximate first order rate constant 
of 5 gmole/hr/g cat. for the 400°F data and the catalyst properties supplied 
by the manufacturer (Appendix A), an effectiveness factor of 83.5% was cal- 
culated meaning that no significant concentration gradients would occur 
within catalyst pellets. 

Examination of a catalyst pellet that had been sliced in two (see 
Figure 4.11) revealed a dark thin outer crust surrounding the inner alumina 
support (white). It is believed that this thin outer layer contained all of 
the ruthenium metal, thus lending support to the high calculated effectiveness 
factor. 


4.4.7 Effect of Hydrogen Sulfide . - Noble metal catalyst for many re- 
actions are notoriously susceptible to poisoning by sulfur containing materials 
particularly hydrogen sulfide. The ruthenium catalyst under study in this 
project was exposed to relatively large concentrations of H 2 S for short periods 
of time and then the ammonia oxidation activity of the catalyst was deter- 
mined. 

In order to obtain a basis for comparison, the unpoisoned catalyst was 
tested in the integral reactor at 400°F with an inlet ammonia concentration 
of 271 ppm and flow rate of 10.5 gmole air/hr. Using a 6.89 g sample of 
catalyst the outlet ammonia concentration was 144 ppm. This point is shown 
in Figure 4.5 by the symbol © . As can be seen from the figure, the con- 
version is in good agreement with the model predicted rate at this condition. 

A 50 cc sample of reagent grade H 2 S was then injected into the feed stream to 
the reactor over a three minute time interval. Since H 2 S is known to be 
oxidized in the presence of many metals, an unknown amount of the material 
could have been oxidized to sulfur or sulfur dioxide before reaching the 
catalyst. The relatively high concentration of H 2 S was used in order to cause 
at least some unoxidized material to reach the catalyst. Based on the inlet 
flow rate the average inlet H 2 S concentration was about 350 ppm. The reactor 
exit was sampled following the injection and the ammonia conversion was found 
to have fallen from 46.9% for unpoisoned catalyst to 30.2%. This point is 
indicated by the symbol © on Figure 4.5. A second 50 cc injection further 
reduced the conversion to 28.4%. This point is indicated by the symbol (3) 
on Figure 4.5. The two injections together caused a reduction in reaction 
rate of approximately 50%. 

Although exposure to hydrogen sulfide does significantly reduce oxida- 
tion activity it does not entirely deactivate the catalyst even in relatively 
large concentrations. 
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4.4.8 Effect of External Mass Transfer . - In any catalytic reaction 
study it is important to ascertain whether or not there are significant 
concentration gradients between the bulk fluid and the catalyst surface. 

If mass transfer to the catalyst surface is a major resistance, the appar- 
ent activation energy which is measured will be lower than the intrinsic 
activation energy when surface adsorption and reaction are the controlling 
resistances. 

Well established techniques are available for evaluating the importance 
of external mass transfer in fixed beds as was discussed previously in 
section 4.2.1. These methods were used directly in calculations involving 
the integral reactor. 

In the case of the stirred-tank reactor, the data were obtained at a 
rotor speed of 550 r.p.m. This can be converted to a catalyst basket linear 
velocity of 242 cm/sec. Because the gas inside the reactor could be swirl- 
ing at an unknown rate in the direction of travel of the catalyst basket, 
the difference in velocity between the gas and the catalyst baskets is un- 
known. It is this velocity difference, which shall be designated the slip 
velocity, which is important in estimating mass transfer coefficients for 
use in predicting external concentration gradients. 

If a slip velocity is assumed, and it is assumed that this slip velo- 
city is analogous to the true gas velocity in the interstices of a packed 
bed containing 30% voids, then the packed bed correlations can be used to 
estimate external concentration gradients in the backmix reactor. The 
following equation was used in the calculation: 


T^(fe) 2 / 3 ■ <»Re>'"- 51 <4J6) 

M K 


where 


By assuming a slip velocity, can be calculated thus making it possible 

to calculate k^, the mass transfer coefficient. Since at steady state the 

rate of mass transfer through the external boundary layer is equal to the 
rate of reaction, the following equation can be used to calculate the con- 
centration gradient: 


r = k' a(P g -P i ) 



a P T (c -c n .) 


(4.17) 
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where Pj = total pressure 

c = bulk concentration 
9 

c^ = concentration at gas-solid interface. 

Knowing the bulk concentration enables one to calculate the interfacial con- 
centration. Knowing the interfacial concentrations and the experimentally 
determined reaction rates, it is possible to obtain kinetic correlations 
which exclude bulk transport effects. This is what was done in the develop- 
ment of Equations (4.14) and (4.15). In obtaining these equations the slip 
velocity was taken to be equal to the basket velocity (242 cm/sec), thus 
allowing values of c^ to be calculated. Values of c i , r N Q , and r^ H were 

then fit by the least squares procedure to Equations (4.14) and (4.15). 

Table 4.3 shows some typical results of the external mass transfer calcu- 
lations which were made using the 242 cm/sec velocity. As can be seen from 
this table the external concentration gradient is negligible at 242°F but 
increases to a moderate level at 400 and 442°F. The slip factor has very 
little influence on the calculated value of either the activation energy 
or the reaction order. The calculated pre-exponential factor decreases 
with a decreased slip factor. These observations indicate that the activa- 
tion energy for the lumped ammonia adsorption and reaction step is truly 
low (5610 cal/gmole), and is not seriously affected by the temperature de- 
pendency of the bulk transport step should the slip factor in actuality be 
less than one. If the slip factor is in actuality lower than the assumed 
value of 1.0, then the implication is that the proposed Equations (4.14) 
and (4.15) would predict rates which are low. Thus, a reactor which was 
operated with absolutely no external mass transfer resistance would show 
ammonia reaction rates and N 2 0 formation rates higher than those predicted 
by Equations (4.14) and (4.15). 

The equations developed using the above mentioned procedure were 
verified with the results from the integral reactor. The Reynolds number 
in the integral reactor were typically about 35 which results in ammonia 
mass transfer coefficients of 1.15 to 1.20 gmole/hr/cm 2 /atm. 

Low intrinsic activation energies for other oxidation reactions have 
been reported in the literature. Johns (ref. 11) in an article on the oxidation 
of trace amounts of CH 4 at 500 to 600° F presents data from which an activa- 
tion energy of 7,250 cal/gmole can be calculated. Satterfield (ref. 9) cites 
research on the oxidation of benzene in which an activation energy of 
2,000 cal/gmole was obtained. In view of these facts, the activation 
energy reported in this work is not unreasonable. 
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4.5 Experimental Data - Platinum Catalyst 

During the early screening of catalyst, supported platinum was found to 
be relatively active for ammonia oxidation (see Figure 4.1). The catalyst 
used and some of its properties are shown in Table 4.1 and Appendix A. In 
order to compare the ammonia oxidation characteristics of this 0.5% platinum 
to the previously studied ruthenium catalyst, 25 experimental runs were made 
using the fixed bed integral reactor. The results of these runs are shown 
in Figure 4.12 along with the calculated results using the ruthenium model. 

As may be seen the ammonia oxidation rate at 400°F is essentially the same 
with platinum as with ruthenium. Humidifying the air to 23% relative humidi- 
ty at 72°F may have resulted in some small increase in the reaction rate. 

The product distribution from platinum is definitely different than 
from ruthenium. Shown in Figure 4.13 is the exit nitrous oxide concentra- 
tion versus inlet ammonia concentration for the integral reactor runs. For 
comparison the ruthenium model predicted nitrous oxide concentration is 
shown. The rate of nitrous oxide production using the platinum catalyst is 
two to three times greater than for ruthenium (see Figure 4.10 for a com- 
parison) . 

These results indicate that although platinum may be more tolerant of 
humidity, its tendency to produce excessive nitrous oxide makes it a infer- 
ior catalyst to ruthenium for spacecabin applications. 


4.6 Experimental Data - Hopcalite Catalyst 

Of the catalyst investigated Hopcalite (see Appendix A and Table 4.1) 
was the only non-noble metal catalyst that showed promise. This catalyst 
is reported by the manufacturer to consist of mixed oxides of maganese and 
copper. Considerable oxidation data has been obtained for the Hopcalite 
catalyst at temperatures between 121 and 352°C with the air flow rate vary- 
ing between 6.3 and 12.6 moles/hr. A summary of the data collected is 
given in Tables 4.4 and 4.5. 

4.6.1 Theoretical Development . - During preliminary evaluation of the 
Hopcalite catalyst isothermal plots of ammonia removal rates as a function 
of average ammonia concentration appeared linear which suggested that the 
reaction is first order with respect to ammonia concentration. 

To determine the validity of this assumption an Arrhenius plot was 
constructed from the data given in Table 4.4. The assumption of first order 
kinetics leads to the rate equation 


r 


Kc 


ave 


F 

W 


(c 1 -c 2 ) 


(4.18) 


37 





Exit Nitrous Oxide Concentration (PPM) 


28 


26 

24 

22 

20 

18 

16 

14 

12 

10 

8 

6 

4 

2 

0 


0.5% Platinum 
T=400°F 

Flow Rate=l0.5 mole/ 
Catalyst Charge=fi,53 







hr 


■ 



▼ 

9 


















































▼ 



■ 

& 







▼ 

■ 

m 




1 





■ 

iM 








■ 






^ 











^ 








/ 






■ 

■ 




0 20 40 60 80 100 120 140 160 180 200 


Inlet Ammonia Concentration (PPM) 

Figure 4.13 - Nitrous Oxide Production Using Integral 


Reactor and Platinum Catalyst 


39 



TABLE 4.4 


EXPERIMENTAL DATA HOPCALITE 


RUN # 
HOPCALITE 

AIR FLOW 
G-MOLE/HR 

TEMP 

°F 

PPM 

IN 

nh 3 

OUT 

PPM N0 2 
OUT 

PPM N 2 0 
OUT 

12/22/70-1 

12.1 

287 

62 

38 

0.2 

20.2 

12/22/70-2 

7.7 

292 

168 

77 

0.2 

11.3 

12/28/70-1 

7.6 

291 

184 

92 

0.2 

12.3 

12/28/70-2 

12.3 

290 

118 

79 

0.2 

4.9 

12/29/70-1 

12.6 

292 

94 

70 

0.2 

2.8 

9/ 1/70-3 

7.4 

300 

252 

126 

0.5 

NA 

9/ 1/70-4 

7.4 

300 

168 

93 

NA 

NA 

9/ 2/70-1 

6.8 

300 

173 

100 

NA 

14.0 

2/26/71-1 

10.4 

305 

151 

88 

0.0 

NA 

3/ 1/71-1 

10.5 

305 

152 

90 

0.0 

9.2 

3/ 8/71-1 

10.5 

305 

176 

117 

0.0 

NA 

3/12/71-1 

10.5 

305 

164 

101 

0.0 

6.3 

3/17/71-1 

10.4 

305 

159 

98 

0.0 

6.5 

3/19/71-1 

10.4 

305 

179 

112 

0.0 

8.5 

3/22/71-1 

10.4 

303 

189 

104 

0.0 

7.9 

3/23/71-1 

10.4 

305 

140 

93 

0.0 

7.3 

3/24/71-1 

10.4 

300 

138 

94 

0.0 

7.6 

3/26/71-1 

10.4 

305 

152 

97 

0.2 

7.3 

3/29/71-1 

10.5 

302 

149 

94 

0.0 

NA 

3/29/71-2 

10.5 

305 

194 

120 

0.0 

NA 

3/31/71-1 

10.4 

302 

112 

70 

0.0 

NA 

3/31/71-2 

10.4 

302 

91 

55 

0.0 

NA 

4/ 7/71-1 

10.4 

304 

156 

98 

0.0 

NA 

4/ 9/71-1 

10.4 

304 

151 

88 

NA 

NA 

10/13/70-1 

10.5 

310 

no 

76 

NA 

1.0 

12/31/70-1 

10.4 

322 

146 

79 

NA 

NA 

1/ 4/71-1 

10.4 

320 

142 

80 

0.0 

NA 

1/ 7/71-1 

10.7 

323 

141 

84 

0.4 

7.1 

1/ 7/71-2 

10.7 

325 

94 

59 

0.5 

5.2 

1/ 8/71-1 

10.7 

324 

90 

56 

0.2 

4.3 

1/ 8/71-2 

10.9 

325 

58 

28 

0.1 

3.4 

1/11/70-1 

9.5 

327 

121 

68 

0.0 

6.0 

1/11/70-2 

9.3 

325 

170 

97 

0.2 

10.2 

9/ 2/70-2 

6.8 

350 

171 

49 

NA 

27.0 

9/ 1/70-1 

10.3 

380 

245 

112 

0.4 

NA 

9/ 1/70-2 

10.3 

380 

118 

54 

NA 

NA 

10/21/70-1 

10.7 

371 

114 

49 

NA 

10.0 

10/22/70-1 

10.7 

371 

61 

40 

NA 

8.0 

10/22/70-2 

10.7 

371 

116 

71 

NA 

NA 

9/ 3/70-1 

6.4 

393 

104 

25 

0.0 

35.0 

10/23/70-1 

10.5 

405 

136 

99 

NA 

20.0 

10/28/70-1 

10.7 

400 

88 

26 

NA 

NA 
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TABLE 4.4 

RUN # 
HOPCALITE 

continued 

AIR FLOW 
G-MOL/HR 

TEMP 

°F 

PPM NH 3 
IN OUT 

PPM N0 2 
OUT 

PPM N 2 0 
OUT 

10/28/70-2 

10.7 

425 

144 

39 

NA 

NA 

10/29/70-1 

10.7 

430 

130 

66 

NA 

NA 

11/ 4/70-1 

10.7 

445 

114 

36 

0.0 

25.0 

9/ 3/70-3 

6.5 

450 

305 

33 

0.0 

109.0 

9/ 3/70-4 

6.5 

450 

169 

19 

NA 

54.0 

8/31/70-1 

6.3 

460 

273 

8 

1 .5 

NA 

9/ 3/70-2 

6.7 

461 

259 

18 

1 .8 

80.5 

11/ 7/70-1 

10.5 

475 

139 

69 

0.2 

33.0 

9/ 3/70-5 

6.4 

500 

218 

21 

5.0 

NA 

11/10/70-1 

10.4 

502 

118 

17 

4.0 

34.9 

11/11/70-1 

10.2 

520 

111 

17 

6.2 

33.0 

11/11/70-2 

7.4 

520 

166 

4 

10.0 

49.2 

11/13/70-1 

7.3 

523 

152 

6 

9.6 

51 .0 

11/13/70-2 

7.3 

523 

166 

19 

9.6 

51 .0 

11/14/70-1 

13.1 

520 

95 

14 

4.8 

NA 

11/18/70-1 

9.0 

515 

229 

42 

4.3 

NA 

11/20/70-1 

7.6 

524 

262 

21 

7.8 

89.9 

9/29/70-1 

12.0 

665 

118 

5 

31.2 

11.0 

4/14/71-1* 

9.4 

300 

58 

50 

NA 

1 .0 

4/19/71-1* 

9.0 

305 

137 

106 

NA 

3.0 
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TABLE 4.5 

EXPERIMENTAL DATA TIME DEPENDENCY OF HOPCALITE REACTION RATE 


RUN # 
HOPCALITE 

AIR FLOW 
G-MOL/HR 

TEMP 

°F 

PPM 

IN 

NH 3 

OUT 

PPM N0 2 
OUT 

PPM N 2 0 
OUT 

2/17/71-1 

10.5 

303 

179 

6 

0.0 

20.3 

2/17/71-2 

10.3 

306 

157 

49 

0.1 

11.4 

2/19/71-1 

10.5 

306 

135 

48 

0.0 

11.3 

2/20/71-1 

10.5 

303 

171 

75 

0.0 

11.3 

2/22/71-1 

10.6 

303 

175 

88 

0.4 

10.1 

2/24/71-1 

10.5 

304 

160 

69 

0.0 

NA 

2/26/71-1 

10.4 

305 

151 

88 

0.0 

NA 

3/ 1/71-1 

10.5 

305 

152 

90 

0.0 

9.2 

3/ 8/71-1 

10.5 

305 

176 

117 

0.0 

NA 

3/12/71-1 

10.5 

305 

164 

101 

0.0 

6.3 

3/17/71-1 

10.4 

305 

159 

98 

0.0 

6.5 

3/19/71-1 

10.4 

305 

179 

112 

0.0 

8.5 

3/22/71-1 

10.4 

303 

189 

104 

0.0 

7.9 

3/23/71-1 

10.4 

305 

140 

93 

0.0 

7.3 

3/24/71-1 

10.4 

300 

138 

94 

0.0 

7.6 

3/26/71-1 

10.4 

305 

152 

97 

0.2 

7.3 

3/29/71-1 

10.5 

302 

149 

94 

0.0 

NA 

3/29/71-2 

10.5 

302 

194 

120 

0.0 

NA 

3/31/71-1 

10.4 

302 

112 

70 

0.0 

NA 

3/31/71-2 

10.4 

302 

91 

55 

0.0 

NA 

4/ 7/71-1 

10.4 

304 

156 

98 

0.0 

NA 

4/ 9/71-1 

10.4 

304 

151 

88 

NA 

NA 


Catalyst Weight 6.604 grams 
*NA not available 
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where r = ammonia removal rate 

c = average ammonia concentration 
ave 3 

K = reaction rate constant 
Ci = ammonia inlet concentration, ppm 
c 2 = ammonia exit concentration, ppm 
F = air flow rate moles /hr 
W = catalyst weight, grams 

If a material balance is made on a fixed bed plug flow reactor, the follow- 
ing equation is obtained 


Fci - Fc 2 = Kc 


(4.19) 


but the exit concentration can be expressed as a function of the inlet concen- 
tration 


dci 

c 2 = c i + air * 


(4.20) 


If this is substituted into Equation (4.19) and simplified, equation (4.21) 
resul ts 


F^.KC . (4.2D 


This can be integrated through the bed from C! to c 2 obtaining: 


In 


Cl KW 

^I = r 


(4.22) 


This can now be equated to Equation (4.18) and the average concentration 
obtained as 


_ ci-c 2 

c ave ~ TrTc77c7 


(4.23) 


Substituting Equation (4.23) for c aye in Equation (4.18) gives 


(4.24) 


jr- = In Ci/c 2 
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The reaction rate is proportional to the concentration of the reactant 
at the surface of the catalyst. This arises since the reaction takes place 
on the surface of the catalyst, but the surface concentration is controlled 
by the rate of mass transfer of the reactant to the particle surface 


r ■ K act C S - k 9 A < c B - C S> 


where K act 

kg 

C B 

C S 

A 


intrinsic rate constant 

mass transfer coefficient 
bulk ammonia concentration 

ammonia concentration on the catalyst surface 

catalyst surface area. 


Equation (4.25)can be solved for the ammonia surface concentration 


(4.25) 


_ kg A c B 

C S ' (K act T kg A) • 


The design equation is 


Substitute for c<. 



M 

P 

d C B 


Pect +k 9 A 1 

F " 

• 

1 Cl 

kg A c B 
L K act +k 9 A _ 


L K act k 9 A J 



Rearrangement gives 


Cl 


_ W 


l _ _ H 

1 c 2 F 


K act kg A 

Tk^ITaT 


(4.26) 


(4.27) 


(4.28) 


(4.29) 
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Equation (4.22) provides an opportunity to calculate K or K obs from the 

kinetic data. This value can be used to calculate the intrinsic rate constant, 

K . , if Equations (4.29) and (4.22) are equated which results in 
act 


obs 


K act k 9 A 


(4.30) 


Equation (4.30) may be rearranged to 


K 


act 


k 9 A K obs 
k 9 A - K obs 


(4.31) 


By using (4.31) the intrinsic rate constant may be calculated from the mass 
transfer coefficient and the observed rate constant which is calculated from 
Equation (4.22). The results of these calculations are shown in Tables 4.6 
and 4.7. 

The Arrhenius equation relates the effect of temperature on the intrinsic 
rate constant and is given by 


K act = A exp(-E/RT) (4.32) 


where A = preexponential constant 
E = activation energy 
R = universal gas constant 
T = temperature °K . 

A plot of In K . vs 1/T will result in a straight line for first order 

kinetic reactions. The slope of this line is equal to -E/R. The Arrehenius 
plot for the data given in Table 4.4 is shown in Figure 4.14. The line in- 
dicated on the plot was obtained from a least squares fit of the data and 
represents an activation energy of 5900 cal/gmole. The value of the preexponen- 
tial constant is 924.4. The first order rate equation can now be expressed as 
a function of both temperature and ammonia concentration 


r = 924.4 exp(-5900/RT) c NH (4.33) 
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46 Figure 4.14 - Hopcalite Arrhenius Plot 




4.6.2 The Effect of Time on Catalyst Activity . - It was frequently ob- 
served that immediately after replacement of the Hopcalite catalyst with a 
new charge, the catalyst activity appeared greater than normal. In order 

to investigate the phenomena a charge of 6.0039 grams of Hopcalite catalyst 
was allowed to operate continuously for 1222 hours at 151°C and 10.5 gmoles 
air/hr. The results of this test are given in Tables 4.5, 4.7, and 4.8. As 
suspected they show a significantly higher activity during the first 200 hours 
of the catalyst utilization. This fact is further depicted in Figure 4.15 
where the observed rate constant is plotted as a function of time. In order 
to avoid this effect all data presented in Table 4.4 were obtained after at 
least 200 hours of operation. 

4.6.3 Performance Data . - Hopcalite catalyst appears to be as effec- 
tive as any of the ammonia oxidation catalyst tested. The rate of ammonia 
oxidation is comparable with both ruthenium and platinum. Figures 4.16, 

4.17, 4.18, and 4.19 show the performance of the Hopcalite catalyst at 142°C, 
151°C, 162°C, and 191°C, respectively. These figures indicate the exit ammonia 
concentration as a function of inlet concentration for a catalyst bed con- 
taining six grams of Hopcalite. Also shown is the predicted results for am- 
monia oxidation over ruthenium catalyst at identical conditions. Because of 
the linearity of the Hopcalite reaction (first order kinetics) the effective- 
ness of the Hopcalite catalyst is slightly better at concentrations greater 
than 80 ppm in the inlet stream. Unfortunately the product distribution for 
the Hopcalite is not as favorable as for the ruthenium. Figures 4.20 and 

4.21 clearly indicate the greater N 2 0 production for the Hopcalite catalyst 
in comparison to the model prediction for the ruthenium catalyst operating at 
identical conditions. 

The distribution of the nitrogen oxides as a function of temperature is 
shown in Figure 4.22. As indicated the conversion of ammonia to nitrogen 
dioxide increases with temperature. Conversely the conversion of ammonia to 
nitrous oxide maximizes at intermediate temperatures. This may lead to er- 
roneous supposition that N 2 0 is a precursor of N0 2 in the reactors. Under the 
comparable conditions no nitrogen dioxide was formed when nitrous oxide was 
substituted for ammonia in the inlet to the reactor. 

4.6.4 The Effect of Mater Vapor on Hopcalite Activity . - To determine 
the effect of water vapor on the activity of Hopcalite catalyst, data was ob- 
tained at 300°F using humidified air. The air was humidified to 23% relative 
humidity at 72°F. The results of this work are presented in Figure 4.23. From 
the data collected on humidified operation it is apparent that water vapor de- 
pressed the ammonia removal rate by approximately 50% for Hopcalite catalyst 

at 300°F. 
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TABLE 4.6 


REACTION RATE CONSTANTS-HOPCALITE 


RUN # 
HOPCALITE 

K obs 

K act 

CO 

o 

X 

0 

j— 

12/22/70-1 

0.976 

1.045 

2.410 

12/22/70-2 

0.959 

1.037 

2.394 

12/28/70-1 

0.841 

0.900 

2.397 

12/28/70-2 

0.799 

0.844 

2.400 

12/29/70-1 

0.618 

0.644 

2.394 

9/ 1/70-3 

0.847 

0.907 

2.368 

9/ 1/70-4 

0.724 

0.768 

2.368 

9/ 2/70-1 

0.621 

0.654 

2.368 

2/26/71-1 

0.945 

1.012 

2.353 

3/ 1/71-1 

0.924 

0.988 

2.353 

3/ 8/71-1 

0.714 

0.752 

2.353 

3/12/71-1 

0.847 

0.900 

2.353 

3/17/71-1 

0.837 

0.889 

2.353 

3/19/71-1 

0.804 

0.852 

2.353 

3/22/71-1 

1.022 

1.101 

2.359 

3/23/71-1 

0.703 

0.739 

2.353 

3/24/71-1 

0.675 

0.709 

2.368 

3/26/71-1 

0.778 

0.823 

2.353 

3/29/71-1 

0.804 

0.851 

2.362 

3/29/71-2 

0.830 

0.882 

2.362 

3/31/71-1 

0.809 

0.858 

2.362 

3/31/71-2 

0.865 

0.920 

2.362 

4/ 7/71-1 

0.808 

0.857 

2.356 

4/ 9/71-1 

0.926 

0.990 

2.356 

10/13/70-1 

0.648 

0.678 

2.338 

12/31/70-1 

1.015 

1.072 

2.302 

1/ 4/71-1 

0.948 

1.015 

2.308 

1/ 7/71-1 

0.890 

0.948 

2.299 

1/ 7/71-2 

0.804 

0.851 

2.293 

1/ 8/71-1 

0.804 

0.851 

2.296 

1/ 8/71-2 

1.263 

1.382 

2.293 

1/11/70-1 

0.888 

0.948 

2.287 

1/11/70-2 

0.845 

0.900 

2.293 

9/ 2/70-2 

1.419 

1.594 

2.222 

9/ 1/70-1 

1.340 

1.473 

2.143 

9/ 1/70-2 

1.351 

1.486 

2.143 

10/21/70-1 

1.504 

1.672 

2.166 

10/22/70-1 

0.775 

0.795 

2.166 

10/22/70-2 

0.874 

0.928 

2.166 

9/ 3/70-1 

1.537 

1.743 

2.110 

10/23/70-1 

0.564 

0.586 

2.081 

10/28/70-1 

2.171 

2.531 

2.093 
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TABLE 4.6 


continued 


RUN # 
HOPCALITE 

K obs 

K act 

1/T °K 

10/28/70-2 

1.925 

2.200 

2.034 

10/29/70-1 

1.207 

1.308 

2.022 

11/ 4/70-1 

2.063 

2.377 

1 .989 

9/ 3/70-3 

2.410 

2.933 

1.978 

9/ 3/70-4 

2.349 

2.843 

1.978 

8/31/70-1 

3.764 

5.222 

1.957 

9/ 3/70-2 

2.985 

3.816 

1.954 

11/ 7/70-1 

1.210 

1.310 

1.925 

9/ 3/70-5 

2.514 

3.072 

1.875 

11/10/70-1 

3.333 

4.215 

1.871 

11/11/70-1 

3.237 

4.060 

1.837 

11/11/70-2 

4.652 

6.855 

1.837 

11/13/70-1 

3.878 

5.298 

1.831 

11/13/70-2 

2.619 

3.198 

1.831 

11/14/70-1 

4.227 

5.587 

1.837 

11/18/70-1 

2.563 

3.079 

1.846 

11/20/70-1 

3.227 

4.140 

1.829 

9/29/70-1 

6.296 

9.701 

1.600 


X 10 3 
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TABLE 4.7 

EXPERIMENTAL DATA TIME DEPENDENCY OF HOPCALITE REACTION RATE 


RUN # 
HOPCALITE 

K obs 

K act 

1/T °l 

2/17/71-1 

5.785 

9.707 

2.359 

2/17/71-2 

2.002 

2.32 

2.350 

2/19/71-1 

1 .811 

2.072 

2.350 

2/20/71-1 

1.436 

1.595 

2.359 

2/22/71-1 

1.221 

1. 34 

2.359 

2/24/71-1 

1.462 

1.628 

2.356 

2/26/71-1 

0.945 

1.012 

2.353 

3/ 1/71-1 

0.924 

0.988 

2.353 

3/ 8/71-1 

0.714 

0.752 

2.353 

3/12/71-1 

0.847 

0.900 

2.353 

3/17/71-1 

0.837 

0.889 

2.353 

3/19/71-1 

0.804 

0.852 

2.353 

3/22/71-1 

1.022 

1.101 

2.359 

3/23/71-1 

0.703 

0.739 

2.353 

3/24/71-1 

0.675 

0.709 

2.368 

3/26/71-1 

0.778 

0.823 

2.353 

3/29/71-1 

0.804 

0.851 

2.362 

3/29/71-2 

0.830 

0.882 

2.362 

3/31/71-1 

0.809 

0.858 

2.362 

3/31/71-2 

0.865 

0.920 

2.362 

4/ 7/71-1 

0.808 

0.857 

2.356 

4/ 9/71-1 

0.926 

0.990 

2.356 


X 10 3 
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TABLE 4.8 


TIME DEPENDENCY OF HOPCALITE CATALYST 
ACTIVITY 


RUN it 
HOPCALITE 

ELAPSED TIME 
HOURS 

nh 3 REMOVAL RATE* 
X 10 5 

K obs 

K act 

2/17/71-1 

0 

30.16 

5.78 

9.71 

2/17/71-2 

8 

18.54 

2.00 

2.33 

2/19/71-1 

48 

15.26 

1.81 

2.07 

2/20/71-1 

67 

16.75 

1.44 

1.60 

2/22/71-1 

119 

15.45 

1.22 

1.33 

2/24/71-1 

166 

15.87 

1.46 

1 .63 

2/26/71-1 

214 

11.02 

0.946 

1.01 

3/ 1/71-1 

286 

10.94 

0.925 

0.99 

3/ 8/71-1 

455 

10.32 

0.715 

0.75 

3/12/71-1 

551 

10.98 

0.847 

0.90 

3/17/71-1 

672 

10.55 

0.837 

0.89 

3/19/71-1 

720 

11.49 

0.804 

0.85 

3/22/71-1 

788 

14.57 

1.02 

1.1 

3/23/71-1 

814 

8.09 

0.704 

0.74 

3/24/71-1 

839 

7.72 

0.676 

0.71 

3/26/71-1 

886 

9.56 

0.779 

0.82 

3/29/71-1 

958 

9.62 

0.804 

0.85 

3/29/71-2 

960 

12.81 

0.831 

0 . 88 

3/31/71-1 

1006 

7.22 

0.810 

0.86 

3/31/71-2 

1008 

6.20 

0.865 

0.92 

4/ 7/71-1 

1176 

10.09 

0.808 

0.86 

4/ 9/71-1 

1222 

10.08 

0.926 

0.99 
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Inlet Ammonia Concentration (PPM) 

Figure 4.20 - N 2 0 Production for Hopcalite Catalyst at 157°C 




58 


(Wdd) uo!4eu}uaouo3 0 Z H 
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Figure 4.21 - N 2 0 Production for Hopcalite Catalyst at 1 62°C 
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Figure 4.23 - Effect of Water on Ammonia Removal Rate 



4.7 Estimated Performance of Catalytic Oxidation System 

In this section information will be presented for estimating the per- 
formance of a fixed bed oxidizer using ruthenium catalyst. This catalyst 
has been shown to be the best one investigated with respect to low nitro- 
gen oxides production. 

The most critical variable in determining the performance of a given 
reactor is the temperature. In order to minimize the formation of nitrous 
oxide, the temperature in the reactor should be as low as is consistent with 
an acceptable ammonia removal rate. This temperature will probably be less 
than 300°F for normal operation. 

The bed geometry does not influence the conversion so long as the 
Reynolds number is kept above about 35. Thus, pressure drop and mechanical 
considerations will determine the shape of the bed. 

The amount of catalyst per mole of inlet air required will depend on 
the inlet ammonia concentration, the desired ammonia outlet concentration, 
the operating temperature, and the presence of poisons in the feed. Figure 
4.24 has been prepared to show graphically the relationship between inlet 
and outlet concentration at various reactor temperatures. This figure can 
be used to calculate the catalyst required to treat any inlet ammonia con- 
centration between about 10 ppm and 100 ppm. For example, at 250°F with 
an inlet ammonia concentration of 50 ppm 2.6 g of catalyst would be required 
per gmole air/hr of feed in order to reduce the outlet concentration to 
10 ppm. At a temperature of 300°F 1.55 g of catalyst would be required to 
reduce a 50 ppm inlet to a 10 ppm exit at a flow rate of 1 gmole air/hr. 

For a larger flow rate, the catalyst charge would be multiplied by the 
flow rate in gmole per hour to get the required catalyst charge. (Note 
again that the particle Reynolds number should be above 35 to keep mass 
transfer resistances small). 

The concentration of nitrous oxide in the reactor exit may be esti- 
mated by the use of Figure 4.25. The nitrous oxide concentration at the 
inlet ammonia concentration is subtracted from the nitrous oxide concen- 
tration at the exit ammonia concentration and the result is the exit con- 
centration of nitrous oxide. For example, at 300°F if the inlet ammonia 
concentration is 50 ppm and the exit ammonia concentration is 10 ppm then 
the exit nitrous oxide concentration is 2.4 ppm. 

In the contaminant control system it is very important that ammonia 
not be subjected to temperatures above about 400°F since large amounts of 
nitrous oxide and/or nitrogen dioxide can result. This means that the feed 
to a catalytic oxidizer used for methane oxidation and operating at 600°F 
to 650°F should be pretreated to remove the ammonia. This pretreatment 
could be a reactor operating at 250°F to 350°F and packed with ruthenium. 
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5.0 EXPERIMENTAL SYSTEMS FOR AMMONIA OXIDATION 


In order to carry out the experimental portions of the ammonia oxida- 
tion studies two reaction systems were utilized. These systems, one a 
backmix (or stirred-tank) reactor and the other a fixed-bed reactor, are 
separately discussed below. The two different systems were selected in 
order to take advantage of the complementing features of having both dif- 
ferential and integral type reaction rate data available for analysis. 


5.1 Backmix Reaction System 

A schematic diagram of the backmix reactor system is shown in Figure 
5.1. Air from a compressed air cylinder passes through a 1/4 in. stain- 
less steel line and into a vertical glass tube containing CaS0 4 dessicant. 
From this point the air flows through a Hastings model LF 3K mass flow- 
meter transducer. The transducer is connected to a meter which provides 
a visual display of the air flow rate to the reactor. After passing through 
the Hastings meter the air passes through a 1/4 in. Whitey needle valve, 
mixes with a small stream of ammonia, and is fed into the reactor. Air 
flow out of the reactor occurs through leaks around the reactor end plates 
and around the agitator shaft. 

Gaseous ammonia is fed to the reactor from a cylinder containing liquid 
ammonia. Ammonia vaporized from the cylinder passes through a pressure re- 
gulator and into a 1/4 in. stainless steel line which contains a 1/4 in. 
Whitey needle valve and connects to a high resistance tube. The resistance 
tube consists of a 3/8 in. O.D. glass tube packed with crushed glass to pro- 
vide high resistance to flow. 

A detailed schematic of the backmix reactor is shown in Figure 5.2. 

The reactor consists of a 4-1/2 in. diameter cylindrical shell, 15 in. in 
length with circular plates at each end. A shaft, which holds the catalyst 
baskets and agitator blades, extends through the centerline of the reactor. 
Thermowells and sample ports are located in both end plates. The reactor 
housing, end plates, agitator shaft, and agitator blades are all constructed 
of stainless steel. The outside of the reactor housing is wrapped with 
electrical resistance heaters which are connected to a variable voltage 
transformer. The outside of the reactor and both end plates are insulated 
to minimize heat losses. The temperature within the reactor is monitored 
with a bare iron-cons tantan thermocouple which is mounted in the end plate 
thermowell and connected to a Honeywell Electronik 15 recorder. With this 
arrangement the temperature can be monitored and controlled anywhere be- 
tween ambient temperature and 900°F. The agitator shaft is driven by a 
variable speed Minarik model SL14 electric motor which is mounted on a 
wooden support which aligns the motor shaft and the agitator shaft. 

The manner in which the catalyst pellets were stacked in the catalyst 
baskets was of concern. The catalyst pellets were cylindrical in shape with 
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Figure 5.1 - Experimental Apparatus 
Backmix Reactor 
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Figure 5.2 - Stirred Tank Reactor Detail 
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a length and diameter of 1/8 in. The first baskets which were made had the 
catalyst stacked one behind the other approximately eight layers deep in the 
direction of rotation. It was suspected for this arrangement that the rear 
pellets were shielded from the gas by the front pellets thereby reducing the 
effectiveness of the catalyst charge. In order to have adequate exposure 
of all the catalyst pellets to the gas, it was necessary to have the catalyst 
spread out so that resistance to flow through the catalyst region was mini- 
mized. This was accomplished using the baskets pictured in Figure 5.3. The 
baskets were made from a flexible steel wire mesh. 

Another area of concern was the degree of mixing within the reactor. 

In order for a backmix type reactor to provide meaningful data it is essential 
that there be no concentration gradients within the reactor. In order to 
check this, runs were made in which samples were taken from both ends of the 
reactor. Good agreement was obtained, at least within the accuracy of the 
analysis, indicating that the mixing was good. As another check on mixing, 
the temperature at both ends of the reactor was monitored. Here again good 
agreement was obtained. 

The fact that the reactor itself was not catalyzing the reaction was 
established. With no catalyst in the reactor, samples were taken from the 
inlet feed line sample port and reactor end plate sample ports. Good agree- 
ment was obtained at temperatures as high as 450°F. At temperatures of 600°F 
and higher, however, some background oxidation of ammonia to nitrogen dioxide 
was noticed. 


5.1.1 Operational Principle . - The backmix reactor is operated in a 
steady state manner in order to take advantage of its differential nature. 

The backmix (or stirred-tank) reactor has a significant experimental advant- 
age over the more conventional gas-solid catalytic reactors such as the fixed 
bed reactor, in that reaction rates as a function of concentration are direct- 
ly calculable. This can be seen by considering the continuous feed backmix 
reactor shown in Figure 5.4. 


F. 

i 



F 


o 


c Ao 


Figure 5.4 - Backmix Reactor Schematic 
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A steady state material balance on component A yields 


F i c Ai * F o c Ao + r A = 0 


(5.1) 


where F 


i 


total molar flow rate in 


F = total molar flow rate out 
o 


u Ai 


inlet concentration of A 


c. = outlet concentration of A 
Ao 

r A = reaction rate of A. 

If the difference in inlet and outlet molar flow rates is zero or negligi- 
ble, then Equation (5.1) can be rewritten 


r A = F i (c Ao -c Ai) 


(5.2) 


If a compound B is being produced by the reaction and no B is present in 
the inlet to the reactor, then the following holds: 

r B * F i c Bo (5 3) 

where c Bq = outlet concentration of B 

r B = production rate of B. 

If the mixing in the reactor is good, it is apparent that the composition 

of the outlet, F , will be identical to that of the reactor contents. Thus 
0 

the reaction occurs at the constant outlet concentrations c^ Q and c B q . 

Relating Equations (5.2) and (5.3) to the specific case of ammonia 
oxidation and nitrogen oxides formation, the following can be written: 


r NH 3 " F i (c NH 3i " C NH 3( 
r N z 0 = F i c N 2 0 


) 


(5.4) 

(5.5) 


r N0 2 = F i c N o 2 

where c MU = inlet ammonia concentration 

NH 3 .j 

c, IU = outlet ammonia concentration 

NH 3q 


(5.6) 
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c N 2 0 = out ^ et N 2 0 concentration 

c NOz = outlet N0 2 concentration 

r NHj = ammonia disappearance rate 
r N 2 0 = n 2 ^ formation rate 

r N0 2 = f° rmat '' on '"ate. 

5.1.2 Experimental Procedure 

Prior to making each run the reactor temperature was adjusted to the 
desired level by manually adjusting the setting of the variable voltage 
transformer. The air flow to the reactor was set to the desired rate by 
adjusting the needle valve downstream of the mass flowmeter. The ammonia 
flow rate was set by adjusting the needle valve between the amnonia cylin- 
der and the resistance tube and observing the pressure on the upstream side 
of the resistance tube. The. purpose of this observation was to get a qual- 
itative idea of what the ammonia concentration at the inlet to the reactor 
will be; for calculational purposes the inlet ammonia concentration was de- 
termined by colorimetric analysis. The agitator was set to the desired 
speed (550 r.p.m.) by adjusting the Minarik motor controller. 

After the temperature and flow rates had been at their steady state 
values for approximately one hour, sampling was begun. The reactor con- 
tents were sampled first. The desired number of samples (usually three 
each) for ammonia, N 2 0, and N0 2 analysis were taken via syringe at the . 
sample port in the reactor end plate, being careful that the rate of sample 
withdrawal did not exceed the air flow rate into the reactor. This was 
important because a rapid sample withdrawal rate could result in room air 
leaking into the reactor thereby diluting its contents. After sampling of 
the reactor contents was completed, the air feed stream to the reactor is 
sampled for NH 3 via a septum in the air feed line downstream of the air- 
ammonia mixing point. 
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5.2 Fixed- Bed Reaction System 

A fixed-bed reactor was constructed in order to provide a system for 
evaluation of ammonia oxidation catalyst. The equipment consists of an iso- 
thermal oven, a feed preparation system, and the reactor proper. 

The oven was constructed of 18 gauge 314 stainless steel and insulated 
with one inch of glass fiber type insulating material. Oven temperature was 
manually controlled by adjustment of the applied voltage to a 500 watt elec- 
trical resistance. The oven temperature was monitored using three thermo- 
couples placed inside of the oven. A long shafted 1/6 horsepower electric 
motor was used to provide air circulation inside of the oven and thus insured 
a uniform temperature distribution. 

A schematic flow diagram of the fixed-bed apparatus is shown in Figure 
5.5. The feed preparation system consisted of an ammonia supply tank connect- 
ed through a high-resistance glass tube to the air feed line. This tube 
provided sufficient resistance to flow to allow very small amounts of ammonia 
to be continuously mixed with the air stream. The air feed was humidified 
prior to the ammonia connection by bubbling through a temperature controlled 
water bath. The humidified and ammonia containing stream then passed by a 
sample port and into the oven. The oven contained a 10 ft. by 1/4 in. stain- 
less steel preheater coil which brought the feed up to reaction temperature. 
The stream then passed upward through the reaction section and out of the 
oven. An exit sample port was located on the gas stream following the re- 
actor. The stream then passed through a wet test meter and was vented. 

The reactor detail is shown in Figure 5.6. The 
of a 1 in. outside diameter type 347 stainless steel 
long. The inside diameter of the tube was 0.765 in. 


reactor was constructed 
tube approximately 6 in. 


71 




Figure 5.5 - Fixed-Bed Reactor System 
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Figure 5.6 - Fixed-Bed Reactor Detail 
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6.0 PERFORMANCE OF AN AMMONIA ADSORBER 


In order to investigate the performance of an ammonia adsorber it is de- 
sirable to have available at least the outlines of a mathematical model of 
the system. This assures that one has a reasonable grasp of the physical pro- 
cesses involved and it reduces the possibility of failure to measure a crucial 
parameter. In addition, these models are very useful in scaling up the system 
and in estimating its performance under conditions not too different from those 
used in testing the model. 


6.1 Modelling Techniques 

When a mathematical model is needed for an engineering system, the first 
approach frequently involves the application of the equations of change. For 
an isothermal ammonia adsorber this reduces to the description of the system 
using an unsteady-state mass balance. In addition to the mass balance it is 
essential to have available a model or relationship for the equilibrium be- 
havior of the system. Since the nature of this behavior is at least as im- 
portant as the differential mass balance, the form of this relationship will 
be discussed first. 

6.1.1 Equilibrium . - What is needed for an ammonia adsorber is an equilib 
rium relationship between concentration or partial pressure of adsorbate in the 
gas phase, and the adsorbate loading on the solid surface. Such a relationship 
when measured and expressed at constant temperature, is called an isotherm and 
a considerable amount of work in the past has been devoted to defining iso- 
therms for various solid-gas systems. Plots of adsorbate bulk concentration 
(or partial pressure) versus equilibrium surface concentration at a particular 
temperature are also called isotherms. 

These isotherms are particularly useful in design because they provide a 
means of estimating the minimum amount of material that would be required to 
adsorb a particular load. In addition, if the breakthrough occurs over a 
narrow time span, they can be used to estimate the time at which the adsorbent 
will break through the bed. 

The Brunauer, Emmett, and Teller (or BET) equation (ref. 12) is, because 
of its easy use and accuracy, probably the most widely used isotherm equation. 
The BET equation can be written as: 


' “m Set 


(Cbet-D p 

"m Set p „ 


( 6 . 1 ) 


where 


C BET 

P 


= dimensionless experimental constant 

= partial pressure of gas 
= vapor pressure of gas 
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Wg. = equilibrium adsorbent loading 

W M = loading, corresponding to a monomolecular coverage of adsor- 
M bent by adsorbate. 

The theory behind the BET equation has several deficiencies. The most 
obvious of these is that the BET model allows an adsorbate molecule to be 
adsorbed on top of an isolated adsorbed molecule. Also, the energy change 
upon adsorption is considered to be equal to that when the vapor is condensed 
to a liquid state in which each molecule has twelve nearest neighbors. De- 
spite these shortcomings, results obtained with the equation have been quite 
good. 


A simplification of the BET equation to the Langmuir equation is possible 
under certain conditions. If P is much smaller than P o (P < 10' 3 atm and P Q > 

10 atm), as is the case for the air-ammonia system, then: 


P 


o 


= Pr 


P 


so that 


if 


and 


then 


P = 


W M C BET 


w e (c bet - 1 ) 


A' 


W M C BET 



B' 




( 6 . 2 ) 


( 6.3 


( 6 . 4 ) 


( 6 . 5 ) 


( 6 . 6 ) 
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which is identical in form to an equation proposed by Langmuir. Since P is 
proportional to C, equation 6.6 can be written as: 


C = 


W E 

A A " W E B B 


(6.7) 


where and B g are experimentally determined constants. 

For the copper sulfate treated sorbeads, as well as the wet and dry 
Hambeads data, an isotherm of this form was used to model the equilibrium 
data. 


Once a suitable form has been determined for the equilibrium relationship, 
it may also be used to gain additional insight into the nature of the adsorp- 
tion that is occurring through the calculation of isosteres. 

6.1.2 Isosteres . - A plot of the natural logarithm of adsorbate partial 
pressure against inverse absolute temperature with the amount adsorbed held 
constant is called an isostere. Because of the difficulty in holding sorbent 
loading constant, isosteres are normally plotted from a fit of isotherm data 
and not directly from experimental results. The procedure for accomplishing 
this will be demonstrated in the results section. 

Isosteric plots are generally straight lines. The slope of these straight 
lines is equal to the negative of the ratio of the latent heat of adsorption 
for the system to the gas constant. An isostere with a large negative slope 
represents a material with a higher heat of adsorption than an isostere with 
a smaller negative slope. 

The Clausius-Clapeyron equation can be written: 


“d" P o> q 

3TT7TT Rg 

where q = latent heat of a phase change, Btu/mol 
Rg = gas constant, Btu/mol - °R 
T' = absolute temperature, °R. 


( 6 . 8 ) 


From equation 6.8, the slope of the natural log of vapor pressure plotted 
against inverse absolute temperature is the negative ratio of the latent heat 
of a phase change to the gas constant. 


Therefore, from an isosteric plot, one can calculate the latent heat of 
adsorption for the adsorbate under study. One can also calculate the latent 
heat of condensation for the adsorbate (care should be taken here to calculate 
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the latent heat of condensation for temperatures at which the vapor pressure 
is near the partial pressures used in the adsorption study). Once the latent 
heat of adsorption and the latent heat of condensation are known, a comparison 
of the two will show whether physical adsorption or chemisorption has occurred. 
That is, if physical adsorption has occurred, the two latent heats will be 
nearly the same. If both chemical adsorption and physical adsorption have 
occurred, the latent heat of adsorption will be greater than the latent heat 
of condensation. 

Although the calculation of isosteres is an important use of the iso- 
therm, the primary function of the adsorption isotherm in this study was in 
conjunction with the differential mass balance and the adsorption rate ex- 
pression. Since this combi nation is a set of relatively intractable equations, 
a numerical procedure was used to calculate the predicted concentration history 
in the bed effluent. 

6.1.3 Gas Phase Material Balance . - The conservation equation can be 
developed for an adsorber by writing an unsteady state gas phase material 
balance for a differential element of the adsorbent bed. For this differen- 
tial element one can write a general balance in the form: 


in - out + generation = accumulation 


Assuming the gas velocity is constant across the cross section of the bed, 
one can write for a differential element: 


in =(C V it R 2 ) 

z 


out =(C V tt R 2 ) 

z+Az 

generation = -ft tt R 2 Az 

AC 

accumulation = tt R 2 Az ^ 

When these terms are substituted into the general balance: 

(C V TT R 2 ) - (C V TT R 2 ) - ft TT R 2 Az = e B TT R 2 Az 

z z+Az 
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rearranging 


(C V) - (C V) 


■^ = ft + 


, AC 
e B AT 


If Az and AT approach zero, the first and third terms of this equation be- 
come the definitions of the partial derivatives so that the material balance 
becomes: 


- IrV * ♦ ‘b If > z • 

If the velocity is also assumed to be independent of axial position in 
the bed, this leads to: 


-v m + 


9C 


9Z 


B 9T 


(6.9) 


6.1.4 Rate of Adsorption . - The term ft in equation 6.9 can be expressed 
in many forms depending on the adsorption model being proposed. However, the 
defining expression for the rate of adsorption is written in terms of the 
adsorbate loading on the solid. The form of the expression used in this work 
is: 


ft = Pfi |f) • (6.10) 

J.n order to be able to use the defining expression for the adsorption 
rate R, it is necessary to determine which of several possible resistances are 
controlling. The first of these to be considered was the resistance to mass 
transfer from the bulk gas phase to the surface of the adsorbate. The applic- 
able form of the general mass transfer equation relative to this situation is: 


- K 
0B 8 


A„ (C-CJ 


( 6 . 11 ) 
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The most difficult term to measure in this expression is the mass trans- 
fer coefficient, K g . However, several good empirical correlations do exist 

which permit reasonable estimates to be made. 

The second resistance which could be controlling concerns the transport 
within the adsorbent particles. If the movement of the ammonia into the pores 
of the adsorbate controls the rate of adsorption, then it would be necessary 
to consider pore diffusion. 

Antonson and Dranoff (ref. 13), using the system ethane-air-molecular 
sieves, have divided the pore diffusion step into macropore and micropore dif- 
fusion (using approximately the same definitions of pore sizes as Dubinin (ref. 

°° 

14), i.e. intermediate pore sizes from radius 15A to 2000A and macropores with 

o 

radius over 2000A). In their work, they concluded that micropore diffusion 
controlled the adsorption rate. 

Descriptive equations of pore diffusion can be developed by writing a 
dynamic gas phase material balance for the pore region. This balance, or 
Fick's law as it is sometimes called, can be written as: 


r 



ac 

-*> 

3r 


3C n 

-*) 

3T 


( 6 . 12 ) 


Equation 6.12 is based on spherical particles and a constant pore size. While 
this is normally not the case, approximations made using these assumptions are 
usually acceptable. Note that here, as in the bulk diffusion case, the ad- 
sorption rate is a function of diffusivity (in equation 6.11, K is a function 
of diffusivity). 9 

Equation 6.12 can be simplified by assuming a pseudo steady state, i.e. 
the term 9Cp/9T may be neglected. With this assumption, equation 6.12 becomes: 


r 



9 a C n 

— -£ ) 
r ar 




(6.13) 


This equation, while not pursued further in this work, is another possible 
expression for the rate of adsorption when pore diffusion is the primary re- 
sistance. 

If neither the movement from the bulk phase nor the diffusion into the 
pores is controlling, the main resistance may be due to surface effects. The 
rate of adsorption on the surface is generally considered in terms of the 


Sf 
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fraction of the surface covered or in terms of adsorption spaces or sites. 

In either case, there is a maximum equilibrium amount of adsorbate (or 
maximum number of filled sites) that can be held on the surface (as a func- 
tion of temperature and local gas phase concentration). It was this maximum 
surface concentration that was discussed by Bruauer, Emmett, and Teller (ref. 12) 
in their paper on equilibrium adsorption. This number of adsorption sites 
is directly related to the sorbent loading for a given adsorbent. It is rea- 
sonable then, to say that the driving force for surface adsorption is a 
function of the difference between equilibrium loading and instantaneous 
loading. In equation form this can be written as: 


ft = YK (W E - W) (6.14) 

When this equation was used in conjunction with the differential mass balance 
and a Langmuir expression for the equilibrium loading, the calculated concen- 
tration at the bed exit computed using a numerical algorithm compared favor- 
ably with the values observed experimentally. 


6.2 Numerical Algorithms 

In order to calculate the concentration of ammonia in the bed effluent, 
a numerical algorithm is needed. 

If all three of the previously discussed resistances, (bulk phase mass 
transfer, pore diffusion, and surface adsorption), were significant in a 
particular process, the descriptive equations would be very complex. In order 
to achieve a manageable set, algorithms were developed which assumed that only 
one of these resistances was significant. Calculations made using the in- 
dividual algorithms were then compared to experimental results to determine 
the value of the models as predictive tools. 

6.2.1 Bulk Phase Mass Transfer Limiting Model . - If bulk phase mass 
transfer limits the adsorption rate, the following equations, which were de- 
veloped previously are applicable: 


V 1£) = e —) + ft 

V 9z ' T e B 3T ; K 


ft- _ 3W 
p B W 


* K g A C (C - C S> 


(6.15) 

(6.16) 
(6.17) 
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Also, since it is assumed in this model that the pore diffusion and surface 
adsorption resistances are negligible and therefore, that equilibrium sur- 
face loading is equal to instantaneous surface loading, it can be written 
that: 


c s ■ fl A - B B u a A A - E b b w e < 6 -' 8 

An inspection of these equations shows that they do not have a simple analy- 
tic solution. However, the equations do lend themselves readily to machine 
computation. Numerical work of this sort requires discrete sized increments 
in time and bed distance and, therefore, double subscripts are required for 
the concentration variables. 

Equation 6.15 can be written in backwards difference form as: 


2F< C iJ 


" S-l.j) = AT {C i,j ' C i,0- 1 } + ^1,j 


(6.19) 


Solving for C. ■ 

* »J 


gives : 



c -V + -i c 

S'-l ,j Az AT ,j-1 


_V 

Az 


+ 


e B 

aT 


- ft 


( 6 . 20 ) 


The defining expression for the adsorption rate can also be written in dif- 
ference form: 


<"i.j - 


( 6 . 21 ) 


For this mass transfer limiting case, the rate of adsorption can also be 
expressed as: 


Vo ' K g A c ’B Vo - c s 


1 ,J 


But using the (Equation 6.18) isotherm equation: 

W i,j 


JS i,j A A- B B W i,j 


( 6 . 22 ) 


(6.23) 
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gives 


W. . - W. 

1 » J AT = K A (C. . 

AT 9 c ' 1 ,J 


w i..i , 

A fl- Vi,] 


(6.24) 


solving for C. . gives 


W. . - W. . W. 

i 3 J i >3""^ j. 1 / a i 9 J 

AT *9 A c A, - BpW. 

1 -J ' K g A c 


(6.25) 


Since 


R i,j ' “i.j-i 1 


(6.26) 


equation 6.20 can be written as: 


e B P B 


C i , J 


= C i-l»J Az + C i,j-1 AT ~ AT ^ M i ,j ~ M i,j-1^ 


_V.fB 
Az AT 


(6.27) 


If the assumption can be made, that C. • ,, C. , and W. . , are known at 

some initial z and T then equations 6.25 and 6.27 can be solved for C. . and 
U 1 »J 

W i ,j* 

A gridwork can be constructed with ordinate of time and abscissa of bed 
length (Figure 6.1). One can also imagine a bulk gas concentration axis per- 
pendicular to the plane of the paper. That is, each point on the time-dis- 
tance grid has associated with it a bulk gas concentration (and an adsorbent 
surface concentration). An examination of this grid shows that at any time, 
the bulk gas concentration at the front of the bed is the inlet concentration 
or: 


C 1 ,j ” C inlet 


for all T > 0 


(6.28) 


! 
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Also, initially, for any position in the bed, the bulk gas concentration and 
the surface concentration (loading) are zero or: 


C. , = W. , =0 for T = 0, 0 < z < Total bed length (6.29) 

1,1 1,1 “ ” 

Therefore, all values of C. . and W. . on the ordinate and abscissa of the 

1 »J 1 *J 

computational grid are defined. Taking the grid position 2-2 (Figure 6.1), 

W. i i. C, ... and C. , . are all known and for this 2-2 position, C. . 

1 > J ” I 1 jJ — ^ 1 — I 9 J 1»J 

and W. . can be calculated using equations 6.25 and 6.27. Obviously if the 
■ »J 

2-2 position is known, the 2-3 position can be solved for W. . and C. .. 

1 » J 1 

Thus a march down the bed at constant time is clearly defined. Once the 
end of the bed is reached, 2 ancl W j 2 are known anc * C-j 3 and W j 3 can be 

calculated. This process can be continued until the entire time-length grid 
has been covered. 

A Fortran IV computer program was written based on the above technique. 
Inputs to the program are the experimental parameters such as inlet ammonia 
concentration, gas flow rate, bed height, and run duration. Also required 
as input are such constants as adsorbent particle size, apparent bulk den- 
sity, surface area per gram of adsorbent and the experimentally determined 
Langmuir isotherm constants. The mass transfer coefficient is calculated 
within the program based on a correlation by Satterfield (ref. 9). Other im- 
portant variables such as gas superficial velocity are also calculated 
within the program. Program outputs include the calculated mass transfer 
coefficient and bulk phase gas concentrations, adsorbent surface concen- 
tration and adsorption rate, all as functions of time and distance down the 
bed. Another program output is a machine plot of bulk phase gas concentra- 
tion versus bed length at constant time for several fractions of the total 
run time. 

Of special interest is the portion of the program in which equations. 

6.25 and 6.27 are solved by trial and error. The Langmuir isotherm equation 
is part of the solution and special care must be exercised so that one always 
operates in the first quadrant of the Langmuir equation. C. • and W. . must 

always be positive in a real system. Due to the asympotic nature of the 
equation, it is easy for the solutions to equations 6.25 and 6.27 to be out- 
side the first quadrant, particularly, since these equations are solved by 
trial and error. Solutions to equations 6.25 and 6.27 outside the first 
quadrant can be convergent but physically meaningless. 

6.2.2 Surface Adsorption Limiting Model .- Several other expressions 
relative to the adsorption rate have been proposed (ref. 15). However, the one 
most applicable to this study is the one introduced previously and rewritten 
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here for convenience: 


ft = YK (W e - W) (6.30) 

Equation 6.30 is similar to equation 6.17 in that both contain a driving 
force in terms of a concentration difference and a coefficient which is in- 
versely proportional to the resistance encountered in that particular step. 
However, the term (W g - W) is proportional to the number of available ad- 
sorption sites on the surface and, therefore, the adsorption rate expression 
for this model assumes that surface adsorption offers the only significant 
resistance to adsorption. This rate expression can be used as was equation 
6.22 above. 

One other step is necessary when using equation 6.30 because W £ is a 

function of the local bulk ammonia concentration. Since in this model the 
surface adsorption resistance is assumed to be much greater than the bulk 
diffusion or pore diffusion resistances, it can be written that: 


C = C s (6.31) 

It has been previously stated that C<. is functionally related to the 
surface concentration by the Langmuir equation. Therefore, it follows that: 

ca A 

W E 1 + CBg ( 6 - 32 ) 

where and B g are Langmuir constants. 

New expressions for C. . and W. . can now be developed based on the 

above discussion but using the same computational sequence as before. The 
new expressions are (in computational order): 


W - C i-1,j A A 


E. . 1 + C. , . B„ 

i,J l-l ,J B 


(6.33) 


W i,j = 


Pg 

YK W F + W . . . 
E i,j AT 


m B 

AT 


(6.34) 


+ YK 
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(6.35) 


1 »J 


c JL+ c -i 

S-1J Az i ,j-l AT 




V . fc B 
Al + AT 


As presented, equations 6.33, 6.34, and 6.35 can be solved in a straight- 
forward manner without the trial and error that was necessary with the bulk 
diffusion model. 


6.3 Experimental Data 

In order to use the numerical algorithms it is necessary to have numeri- 
cal values for the parameters contained within these models. Since the 
equilibrium relationship plays a key role in both of these algorithms, the 
experimental method used to evaluate the two constants contained in the 
Langmuir form of this expression will be discussed first. 

6.3.1 Bed Capacity Calculations . - A graphical integration technique 
was used to calculate adsorbent capacities. If the plot of exit ammonia 
concentration versus time is graphically integrated and the resulting nunfoer 
divided by total run time, the time averaged exit ammonia concentration is 
obtained. Since the inlet ammonia concentration is held constant, the dif- 
ference between the inlet and average exit concentrations is the average 
ammonia concentration adsorbed. This concentration in mols per ft 3 can be 
multiplied by the total run time and the flow rate to give the moles of am- 
monia removed during a run. Since the amount of material in the bed is 
known, the bed capacity in grams of ammonia per gram of bed may be calculated. 

Bed saturation runs were made at temperatures of 72°F and 100°F. Air 
flow rates ranged from 0.164 cfm to 0.202 cfm while ammonia inlet concen- 
trations went from 90 ppm to 520 ppm and bed length ran from two to eight cm. 
Fresh adsorbent was used in all runs except Run 24R. 

The sorbent for Run 24R was the same sorbent used in Run 24. The 
material from Run 24 was regenerated in an oven at 383 ± 5°F for five days. 

A plot of exit ammonia concentration against time for Run 24 and 24R is 
shown in Figure 6.2. 

The gravimetric capacity listed in Table 6.1 is the difference of the 
bed weight after and before a run divided by bed weight before a run. From 
Table 6.1, it can be seen that calculated bed capacities are uniformly higher 
than the gravimetric capacities. This, coupled with an observed color inter- 
face which passed through the bed during the run and the fact that no com- 
pounds other than air and ammonia were detected (no attempt was made to 
analyze for water) in the bed effluent, leads one to consider that some chemi- 
sorption was taking place in the bed. A specific chemical reaction which 
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Figure 6.2 - Experimental Points 
for Runs 24 and 24R 
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Pressure drop across bed 0.2 in Hg guage. 

Run 23 scrapped because of experimental error. 
CSTS - Copper Sulfate Treated Sorbeads 
DHMB - Dry Hambeads 
WHMB - Wet Hambeads 


would give such results is the replacement 

CuS0 4 • 5H 2 0 + 4NH 3 t CuS0 4 • 4NH 3 + 5H 2 0 

The color change could be accounted for by the fact that the characteristic 
color of CuS0 4 • 4NH 3 is dark blue which was the color that the bed turned. 

Adsorbent from Run 24 was regenerated at 385°F and resaturated as Run 
24R. At these rather mild conditions the adsorbent was regenerated to ap- 
proximately 80 percent of the fresh adsorbent's capacity based on calculated 
bed capacities. At more severe conditions, perhaps 500°F, the adsorbent 
would likely be regenerated to an even higher fraction of the original capa- 
city. The regeneration, color change, and gravimetric results leads one to 
suspect that both chemical and physical adsorption were taking place in the 
bed. 


6.3.2 Determination of the Languir Constants . - When there is no change 
in concentration of the ammonia stream as it passes through the bed then the 
adsorbent within the bed is in equilibrium with that stream. Therefore, the 
calculated values for the bed capacities can be used in conjunction with a 
least squares analysis to find the best fit of the data to an equilibrium re- 
lationship of the Langmuir form. 

Figures 6.3 through 6.5 illustrate the agreement between the fitted ad- 
sorption isotherms and the experimental results. As can be seen this two 
parameter model is adequate to describe the data over the entire experimental 
range. The form of the Lanpiuir isotherm equation is such that curves drawn 
from it pass through the origin (as would be expected on theoretical grounds 
for all isotherms). 

6.3.3 Isostere Analysis . - Isosteres were plotted (Figure 6.6) using 
the copper sulfate treated sorbeads at three adsorbent loadings. The data 
for the plots were obtained from the fitted isotherm equations and not from 
the experimental points. (The loadings used were, however, within the range 
of the experimental sorbent loadings). The data for the plots were generated 
using equation 6.7. This equation is valid for both 72°F and 100°F (with the 
proper temperature dependent constants). Therefore a given sorbent loading 
and the appropriate constants can be substituted into the equation and the 
free stream concentration at both temperatures can be calculated. 

In Figure 6.7, ammonia vapor pressure is plotted against inverse abso- 
lute temperature on semi -log coordinates at vapor pressure of from five to 
60 mm Hg. The data for this plot is taken from thermodynamic charts for 
ammonia. 

Since the latent heat of adsorption is directly proportional to the 
negative slope of the isosteric plot, this latent heat of adsorption can 
easily be calculated from Figure 6.6. Also, since the latent heat of con- 
densation is directly proportional to the negative slope of the vapor pres- 
sure plot, this latent heat can be calculated from Figure 6.7. These 


89 









Figure 6.5 - Least Squares Fit of Wet Hambeads 
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Figure 6.7 - Plot of Ammonia Vapor Pressure 


vs Inverse Temperature 



calculations indicate an average latent heat of adsorption of about 40,000 
Btu/mol and a latent heat of condensation of about 11 ,000 Btu/mol . This 
means that, because of the higher heat of adsorption, chemical reaction is 
probably taking place during the adsorption and is, therefore, contributing 
to the adsorption. This supports the conclusion that the sorption capacity 
of the bed may not be determined gravimetrical ly. 

6.3.4 Physical Properties of the Adsorbents . - In addition to the 
equilibrium relationship, the density and the void fraction of the sorbent 
in the bed must be known before the algorithms can be used. When the mass 
transfer rate in the gas phase is assumed to be limiting, the external area 
of the adsorbent must also be measured experimentally for substitution into 
the algorithm. 

The void fraction of the bed was determined using benzene as the dis- 
placing liquid. A sample of adsorbent was placed in a graduated cylinder 
and the level volume recorded. The benzene was added (Benzene was used 
because it does not adsorb on silica gel strongly enough to cause thermal 
cracking of the gel.) to some level in the cylinder (above the adsorbent). 
The difference in volume between the adsorbent and benzene was then sub- 
tracted from the total amount of benzene added. This adjusted volume of 
benzene added, divided by the volume of adsorbent in the cylinder, repre- 
sented the total void fraction of the adsorbent. The experiment was re- 
peated before the benzene in the adsorbent's internal pores has evaporated, 
and, therefore, the void fraction calculated for the second experiment was 
the external void fraction. This procedure was repeated several times and 
an average void fraction was used in subsequent calculations. 

Some physical data concerning the surface area of the adsorbent was 
mpacurpH for rnnDer sulfate treated silica gal (CSTS) but not for Hambeads. 
The Hambeads were received at a later point in the study when preliminary 
testing of the algorithms had indicated that these parameters did not need 
to be measured. Table 6.2 contains a summary of this data. 


6.4 Testing of the Numerical Algorithms 

One of the most important components of a model for the system was 
developed when the adsorption isotherms were adequately described using 
an equation of the Langmuir form. Once this relationship had been modelled, 
it could be used in a testing program of the nuemrical algorithms. 

6.4.1 Bulk Phase Mass Transfer Model . - The first model tested was the 
bulk phase mass transfer controlling model. In this model it was assumed 
that pore diffusion and surface adsorption resistances to adsorption are in- 
significant relative to the overall adsorption rate. This model can be re- 
duced to two equations and two unknows as shown by equations 6.25 and 6.27. 
Therefore, this model has a straightforward solution. 

A plot of exit ammonia concentration versus time as predicted by the 
mass transfer model is given in Figure 6.8. Superimposed on this plot is a 
set of typical experimental data from Run 21. From this figure, it is seen 
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TABLE 6.2 


Total void fraction 
(measured) 

External void fraction 
(measured 

Density 

(measured) 

Particles per gram 
(measured) 

External area per gram 
(calculated) 

Particle "diameter" 
(measured) 


ADSORBENT PROPERTIES 

CSTS 

0.61 ± 0.03 

0.325 ± 0.02 

51. ± 0.2 lbs /ft 3 
m 


Hambeads 

0.52 ± 0.06 

0.38 ± 0.01 

52. ±1.0 lbs ffl /ft 3 


313. ± 10 

0.0717 ± 0.0002 ft 2 /g 
0.0724 ± 0.004 in 
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that the mass transfer model predicts a much faster saturation of the ad- 
sorbent than actually occurs. This is not surprising when one considers 
the computed magnitude of the mass transfer coefficient. A mass transfer 
coefficient of 1000 - 1200 ft/hr is considered quite large by most sources 
and would lead to an early saturation of the bed. 

Since the computed mass transfer coefficient is subject to some error, 
the model was also run for values of 120 percent and 80 percent of the com- 
puted mass transfer coefficient. Differences in the three computer runs 
were insignificantly small and therefore the breakthrough curves predicted 
by the runs with high and low mass transfer coefficients were not shown in 
the figure. 

Useful information can, however, be gleaned from this model. First, 
both the size of the mass transfer coefficient and the shape of the pre- 
dicted breakthrough curve indicate that the resistance to bulk phase mass 
transfer is small and is probably not a significant factor in the overall 
adsorption process. Also, the basic modelling procedure is sound and con- 
vergent. This was determined by running the program several times until 
the results indicated that the model gave consistent answers which are 
neither time nor length dependent. Time and length increments were varied 
over a range of one half to twice the basic time and length increments with 
negligible change in the final results. 

The trial and error nature of the solution to the mass transfer equations 
(Equations 6.25 and 6.27) necessitates using computer search techniques not 
generally needed in adsorption calculations. Since this procedure was 
required in only the mass transfer model, it was incorporated in a subpro- 
gram in the general computer program. 

6.4.2 Surface Adsorption Model . - Experimental runs were made with two 
different types of adsorbent in the bed, Hambeads and copper sulfate treated 
sorbeads (CSTS). The humidity level in the feed was zero for all of the 
CSTS runs, but three of the Hambead runs were made with the relative humidity 
of the feed at approximately 50%. 

Figures 6.9 through 6.12 are plots of the calculated ammonia exit con- 
centration for the surface adsorption controlling case versus time for all 
of the runs made with copper sulfate treated sorbeads except Runs 24R and 36R. 
Superimposed on the plots are the experimental data for those runs. The agree- 
ment is quite good on all runs 1 except Run 25. This agreement between cal- 
culated and experimental points suggests that, over the range studied, surface 
adsorption, by the proposed mechanism, is the controlling resistance to ad- 
sorption. 


J The constant YK of equation 6.30 was determined by trial and error. In 
this trial and error, the sum of differences between calculated and experi- 
mental points (at constant time) squared was minimized. The YK shown in the 
figures was selected when YK's 25 percent higher and 25 percent lower gave 
higher sums of squares. A five percent change in YK produced a negligible 
change in the shape of the breakthrough curve. 
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From Figure 6.11, it can be seen that the observed adsorption rate for 
Run 25 is higher than the predicted rate. Aside from this, the predicted 
breakthrough curves have the same basis shape as those observed experimentally. 

For all of the runs conducted at 72°F, a value of 20 lb/ft 3 hr for the 
constant YK produced an acceptable correspondence between the predicted and 
observed bed effluent concentration. At 100°F, the preferred value of YK 
changed to 80 lb/ft 3 hr. 

Runs 29 and 31 were made at 100°F with Hambeads in the adsorber and no 
moisture in the gas stream fed to the bed. The best least squares value for 
YK in these two cases was 20 and 29, respectively, with 25 being a reasonable 
value to use in both cases. A comparison between the predicted and observed 
concentrations at the bed exit are given in Figure 6.13. 

Only one run, number 35, was made at 72°F with a dry ammonia feed and 
Hambeads in the bed. There is considerable scatter in the experimental data 
which leads to an unusually high least squares value of 54 for YK. However, 
as Figure 6.14 indicates a YK value of 25 still provides a reasonable fit of 
the data. 

Three runs were made at 72°F when the relative humidity of the feed to 
the bed of Hambeads was approximately 50%. These data and fitted curves are 
shown in Figures 6.13, 6.14, and 6.15. In one of these runs, number 37, there 
is a considerable amount of scatter in the experimental data which leads to an 
unusually high least squares value of 81. However, for the other two runs, 
the least squares value of YK is the same at 25. Even for Run 37, a YK value 
of 25 does a creditable job of describing the data. This is shown in Figure 
6.15 where the experimental and computed profiles are compared. 

All of this data suggests that a YK value of 20 to 25 would be a good 
design value when considering adsorbent beds of Hambeads or copper sulfate 
treated silica gel at a temperature of 72°F. This value is probably conser- 
vative when a moisture level of approximately 50% relative humidity is added 
to the stream. 


6.5 Estimated Performance of an Adsorption System 

The performance of any adsorption system will be strongly influenced by 
the nature of the adsorption isotherm for the material in the bed. In Figure 
6.17 a comparison is made between the fitted euqations for all of the iso- 
therms used in this study. As the data in the figure indicates, the copper 
sulfate treated sorbeads had the highest capacity at inlet gas concentrations 
greater than 30 parts per million. No speculations are made concerning the 
reason for this difference in capacity since no information concerning the 
preparation of Hambeads was available. The preparation of CSTS is discussed 
in section 7.0. 
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Figure 6.14 - Comparison of Computed and Observed 
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Figure 6.17 - Comparison Between the Fitted Isotherms 



An estimate of the minimum amount of adsorbent that would be required 
to adsorb a particular load can be obtained using these curves. As an ex- 
ample, assume that a one pound per hour stream which contains 100 parts 
per million ammonia is to be passed through a bed of CSTS at 72°F. For 
one days operation: 


1 lb # 24 hr 
hr * day 


100 X 10 -6 lb NH 3 
lb air 


2400 X 10 -6 lb NH 3 
day 


The capacity of CSTS at 100 ppm and 72°F is 0.0161 lb NH 3 /l b adsorbent. 
Multiplying these two numbers indicates that: 


2400 X 1CT 6 lb NH 3 
day 


lb adsorbent 
0.0161 lb NH 3 


0.15 1b adsorbent 
day 


At 100°F, the capacity of CSTS drops to 0.005 lb NH 3 /l b adsorbent 
and therefore 0.48 pounds of adsorbent would be required as a minimum. 
Using a density of 51 lbs m /ft 3 would allow the volume of these beds to be 
estimated. T 


me date ccr.ccrr 0 ^ >*h the reaeneration of CSTS indicate that, for 
one run at least, 80 percent of the fresh adsorbents capacity woo legair.cd 
upon regeneration at 383°F. This suggests the possibility of having two 
beds with a swing cycle or of having one large bed which is removed from 
service occasionally and regenerated. No attempt has been made to optimize 
the regeneration time or temperature. However exposure to the environment 
in outer space would probably result in drastic reductions in both time 
and temperature. 

None of the adsorbents have the ideal shape for their breakthrough 
curves. If a perfect adsorbent could be ordered for this service, one 
would undoubtedly specify that the ammonia concentration in the effluent 
remain near zero until the entire bed was saturated. A review of the break- 
through curves for CSTS and Hambeads reveals that for the small beds used 
in this study, the effluent attains 50% of the influent concentration when 
approximately 20-25% of the time required to reach steady-state has passed. 
Neither sorbent appears to possess any significant advantage over the other 
relative to this desirable feature. 

The shape of the breakthrough curve can be estimated using the surface 
adsorption limiting algorithm once an estimate for YK is available. As 
stated previously, a value of 20-25 would be a good design value when using 
Hambeads . 
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More data is available for copper sulfate treated sorbeads and so the 
variation of YK with temperature can be estimated. If this rate constant 
is assumed to follow an Arrhenius form in a manner similar to kinetic con- 
stants, an equation of the form: 


YK = H' 


-(E'/RT 1 ) 
e y 


would be in order. The values of YK = 20 at 72°F and YK - 80 at 100°F can 
be used to obtain the two constant H' and E ' : 


H' = 2.2 X 10 13 lb/ft 3 hr 


E' = 29,000 Btu/lb mole 


Since data were collected at only two temperatures, the validity of 
these constants cannot be checked. However Daniels and Alberty (ref. 116) have 
considered similar reactions and obtained values of H 1 and E' of the same 
order of magnitude. 

Due to the reduced capacity of the bed at the higher temperatures 
there is little incentive to raise the temperature above 72°F. This equa- 
tion would be useful however in making estimates of minor changes around 
72°F. 


In summary, the isotherms provide a lower limit on the amount of 
material required while the small amount of regeneration data suggests the 
possibility of a significant retention of capacity upon regeneration. 
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7.0 EXPERIMENTAL SYSTEM FOR AMMONIA ADSORPTION 


In order to measure the capacities of the adsorbents it was necessary 
to construct an experimental apparatus with which the concentration, flow 
rate, and temperature of the feed stream to the bed could be held constant. 
The effluent concentration was measured using the techniques discussed in 
section 8.0. 


7.1 Description of the Adsorption Apparatus 

The apparatus used in this study was simple in concept and design. The 
essential elements were a metered air supply, an air-ammonia mixing section, 
a fixed support for the adsorption bed and a thermostatic bath surrounding 
the adsorber. 

Specifically (Figure 7.1), compressed air at 80 psig from the university 
air system was cleaned by passing it over a four to six mesh silica gel bed 
contained in a carbon steel holder. This holder had an inside diameter of 
approximately two inches and was about four inches long. Air saturated with 
water vapor at 80 psig is relatively dry (2200 ppm water) at the atmospheric 
pressure of the study so the main function of the silica gel was to remove 
dust, rust, and entrained oil. (Bimonthly changes of silica gel were made. 

At no time was more than the upstream half of the bed soiled). Air flow was 
maintained via a small air regulator upstream of the silica gel bed. 

Air flow (0.16 - 0.20 ± 0.05 cfm) was measured by a calibrated mass 
flowmeter (Hastings Mass Flowmeter and Mass Flow Transducer by Hastings 
Raydist, Inc., Hampton, Virginia). Tiii^ ;r.ctcr ccjld Kq and adjusted 

for a zero flow easily, and the meter never required more than a ten percent 
adjustment due to instrument drift. 

Air and ammonia were mixe4 in a tee (all tubing and fittings were 1/4 
inch 316 stainless steel by Swagelok unless otherwise specified) approximate- 
ly 13 linear feet from the adsorber. Prior to mixing, the ammonia was stored 
under pressure in an aluminum tank about two feet high with an inside diameter 
of nine inches. The transfer lines from this tank to the mixing tee were of 
glass. Teflon, and stainless steel. A 1/4 inch three way valve was provided 
just downstream of the mixing tee. 

The fixed support for the adsorber was made of 21.5 - 22.0 mu inside 
diameter glass tubing, an aluminum adsorber support plate and Teflon packing 
under the support plate. An excellent glass on steel seal was provided by 
nominal one inch Teflon Swagelok ferrules at each end of the glass tube. 

Given a few minutes to "set", the seal showed no leaks at the pressure en- 
countered in this study. 

Thermostatic conditions were maintained in the adsorber by a constant 
temperature bath (Sargent Thermonitor Model ST by E. H. Sargent and Company 
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and Heater and Circulator for Thermostatic Baths by Sargent-Wel ch Scientific 
Company). The bath was used only for the 100°F adsorption runs and the tem- 
perature control here was excellent with no noticeable variations. A coil 
of tubing was placed in the bath and upstream of the adsorber to insure that 
the bath temperature was the same as the air-ammonia temperature entering 
the adsorber. 

Valves, pressure indicators (PI), thermometers (TI), and sample ports 
were located strategically throughout the process in order to maintain steady 
flows, and to monitor the temperature and ammonia concentrations into and out 
of the adsorber. 

All materials used in this study were tested and found to be chemically 
inert to parts per million concentrations of ammonia. This was determined 
by passing a known parts per million concentration of ammonia through the 
empty system and ascertaining that the inlet and exit ammonia concentrations 
were the same. 


7.2 System Flows 

Compressed air from the air header flowed first through the regulator, 
through the silica gel bed, (Figure 7.1), then through the mass flowmeter, 
and finally to the air-ammonia mixing tee. From there the air-ammonia mix- 
ture went either to vent or through the adsorber bed depending on the position 
of the three way valve. This procedure allowed the concentration of the air- 
ammonia mixtures to be adjusted before starting an adsorption run. 


f.5 txper Mlieil Ld I riutcuui c 

For at least 48 hours prior to making a run, the copper sulfate impreg- 
nated silica gel was preconditioned by putting it in the adsorber and flowing 
process air over it. This was done so that the silica gel was properly 
acclimated to the air and so that any observed changes during an adsorption 
run were solely due to the presence of the ammonia. Weight changes in the 
silica gel were under 0.02 grams for this preconditioning step. 

Just prior to starting a run (after preconditioning) the air flow was 
diverted to the vent and the ammonia holding tank was pressurized using a 
lecture bottle of liquid ammonia. After the concentration of ammonia in the 
air-ammonia mixture had reached steady state in the vent line, this vent was 
sampled to determine the particular concentration. If the ammonia concen- 
tration was not acceptable for the particular run, ammonia vapor was bled from 
or added to the holding tank (thereby lowering or raising the holding tank's 
pressure). If the ammonia concentration was acceptable, the flow was directed 
into the adsorber, the timer was started, and a sample of the effluent from 
the adsorber was taken. 

Samples of adsorber exit and inlet gas flows were taken (in the ratio of 
four or five to one) on a rapid basis at the start of a run. As the bed 
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became saturated with ammonia, the rate of change of exit concentration 
with time decreased and samples were taken less frequently. 

System variables, such as air flow and bath temperature, were monitored 
and recorded approximately every 60 minutes. A given run was considered 
over when the exit concentration approached within ten parts per million of 
the inlet as this was considered to be within the accuracy of the analysis. 

After a run was completed, the adsorbent bed was removed, weighed, and 
stored. The empty system was then flushed with pure air for at least one 
hour to desorb ammonia from the system and thereby to start the next run with 
an ammonia free system. 

Samples at the bed entrance and at the exit from the adsorber were 
analyzed for ammonia using Nessler's Reagent. Samples were removed by gas 
syringe through silicone rubber septa. The details of this test are given 
in section 8.0. 

A periodic check of the entire gas stream was made using a Varian 
Aerograph gas chromatograph equipped with a column packed with approximately 
nine feet of Porapak Q packing. This chromatographic check was made to de- 
termine if any oxides of nitrogen were being formed in the adsorber bed. 


7.4 Preparation of Copper Sulfate Treated Silica Gel 

Silica gel in the form of Mobil Oil Company's Sorbead R was treated with 
an aqueous solution of copper sulfate, allowed to come to equilibrium with 
the solution, and then dried in an oven. Specifically, the procedure was as 
follows: 


100 grams of CuS0 4 • 5H 2 0 (from batch 34385 of J. T. Baker 
Company) were dissolved in 200 milliliters of distilled 
water to give an unsaturated solution of copper sulfate in 
water (some heating required). About 200 milliliters of 
Mobil Sorbead R was added to the solution. About four 
hours, with stirring, was allowed for the system to come 
to equilibrium. The Sorbeads -copper sulfate solution was 
then vacuum filtered and dried in an oven at 350-450°F. 

Finally the beads were crushed to the 14 to 16 mesh size 
used in the study. 

Untreated Mobil Sorbead R is the characteristic off-white of silica gel. 

The treated beads took on a light green hue. 
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8.0 


CHEMICAL ANALYSIS METHODS 


This study of the removal of ammonia required chemical analysis techni- 
ques for the quantitative measurement of ammonia, nitrogen dioxide, and ni- 
trous oxide in the range from about 5 ppm to 500 ppm in a background of 
humidifed air. In order to be most useful the methods needed to be accu- 
rate, simple, rapid, harmless, and inexpensive. The gas sample size avail- 
able for analysis was generally not more than about 100 cc at ambient tem- 
perature and pressure. After considerable investigation, the analysis 
methods chosen were colorimetric in the cases of ammonia and nitrogen di- 
oxide, and chromatographic for nitrous oxide. 


8.1 Ammonia Analysis 

A relatively extensive effort was made in the early stages in this 
investigation to develop a reliable analysis for ammonia using a helium 
ionization gas chromatograph (Varian model 1532-B). This work is summa- 
rized in the following section. The ammonia analysis method finally 
chosen was a colorimetric technique using Nessler's reagent. A descrip- 
tion of the method is given in Section 8.1.2. 

8.1.1 Chromatographic Ammonia Analysis . - Many different schemes 
for the chromatographic analysis of ammonia were tried during the course 
of the project. Ammonia peaks under all conditions had extremely long 
retention times and tailed very badly. Summarized below are the columns 
and conditions that were investigated. 


Support 

Dimensions 

Temp. Range 

Column Const. 

Porapak Q 

1/8!L x 20’ 

Ambient to 
150° 

SS 

Porapak Q 

1/8" x 10' 

Ambient to 
150° 

SS 

Porapak R 

1/8" x 15' 

Ambient to 
150° 

SS 

Chromo 104 

1/8" x 15' 

Ambient to 
150° 

Teflon Lined 
A1 umi num 

Porapak Q 

1/8" x 12' 

Ambient to 
150° 

Teflon Tubing 
Encased in SS 
Shield 
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The shape of the ammonia peak is indicative of physical adsorption of 
ammonia on the column and chromatographic hardware. This is further 
strengthened by the fact that a five foot length of unpacked tubing for a 
column gave identically shaped peaks. In order to minimize the adsorption 
of ammonia on the chromatograph hardware approximately 50% of the tubing 
and fittings which are outside the oven were wound with electrical resis- 
tance wire and heated to 120°C. This also had little effect on the ammonia 
peak shape. 

Effort was also directed toward using a very short column. This 
would allow everything but ammonia to pass through the column quickly as 
well as minimize the column surface area. This effort did not signifi- 
cantly improve the ammonia peak shape. Further attempts to improve the 
ammonia peak shape were to inject directly into the column and to use 
Teflon coated tubing. Neither of these methods have proved satisfactory. 

8.1.2 Colorimetric Ammonia Analysis and Calibration . - Gas samples 
in this investigation were analyzed for ammonia colorimetrically using 
Nessler's reagent in a technique similar to that described in Standard 
Methods (ref. 17). The method consists of contacting water containing ab- 
sorbed ammonia with Nessler's reagent and reading the transmittance of 
the resulting yellow-brown solution against an appropriate blank on a 
colorimeter. The intensity of the color change which occurs is propor- 
tional to the ammonia concentration. 

The step by step analytical procedure is as follows: 

1. 50 cc polypropylene syringe is loaded with 2 cc of distilled 
and deionized water. 

2. The water is contacted with sample gas by drawing 48 cc of 
gas into the syringe, being careful to let the gas bubble 
through the water. The syringe is shaken for approximately 
one minute to promote ammonia absorption. 

3. The contents of the syringe are transferred to a cuvette and 
1 cc of Betz Nessler reagent No. 246 which has been diluted 
5 to 1 with distilled water is added. 

4. After allowing approximately ten minutes for color development, 
the sample is read on a Bausch and Lomb Spectronic 20 colori- 
meter at 420 my against a blank containing 2 cc of distilled 
water and 1 cc of diluted Nessler reagent. The sample cuvette 
and the blank cuvette are chosen such that they yield the same 
transmittance values when filled with distilled water. 

Approximately midway through the experimentation it was found advanta- 
geous to modify the syringe sampling technique by loading the syringe with 
3 cc of water instead of two and drawing two 47 cc gas volumes instead of one 
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48 cc volume. For this case the blank consisted of 1 cc of diluted Nessler 
reagent and 2 cc of distilled water. The new procedure had the effect of 
decreasing the deviation in percent transmittance readings between samples. 
With modified sampling procedure a standard deviation of 3.4 ppm was cal- 
culated at 121 ppm versus 5.1 ppm at 135 ppm for the original procedure. 

Since there was no interaction between precision and concentration, these 
numbers seem to indicate an increased precision for the modified procedure. 
Three determinations were made of each gas stream sampled. 

In order to have a quantitative analysis it was necessary to obtain a 
calibration curve which relates ammonia concentration to percent transmit- 
tance. This was done by preparing known concentrations of ammonia in air 
and establishing the percent transmittance associated with each concentra- 
tion. Known ammonia concentrations were prepared by passing measured amounts 
of air over a liquid ammonia filled Teflon permeation tube in the apparatus 
diagrammed in Figure 8.1. The permeation tube, when held at constant tem- 
perature, lost weight linearly with time and thus allowed one to know the 
exact concentration of ammonia when the air flow rate was known. The per- 
meation tube weight loss was determined by periodic removal from the ap- 
paratus and weighing. 

The apparatus in Figure 8.1 was operated by passing air from a com- 
pressed air cylinder through a coil in a constant temperature bath, over 
the ammonia filled permeation tube and through a calibrated wet test meter. 
The permeation tube was held at constant temperature in a jacketed glass 
housing. Water from a constant temperature bath was circulated through the 
housing via a positive displacement pump. 

cwm, the ahnve calibration procedure it was possible to prepare the 
calibration curves of Figure 8.2. ine upper line uT riyu.c C.2 \:zz 
using 2 cc of water in the syringe and one 48 cc gas draw, while the lower 
line was obtained with 3 cc of water in the syringe and two 47 cc gas draws. 
The linearity of both lines shows that Beer's law is obeyed in the concen- 
tration range investigated. 

The only interference problem encountered with the Nessler reagent 
was that aldehydes cause the solution to be clouded and give false reading. 
Careful attention to the purity of the water used in preparing reagents 
and avoidance of aldehydes around the analysis area prevented interference. 


8.2 Nitrogen Dioxide Analysis 

The analysis of N0 2 was performed in a manner similar to the ammonia 
analysis. Sample analysis was performed colorimetrically using the modi- 
fied Saltzman reagent (ref. 18) and the syringe method of Meador and Bethea 
(ref. 19). 

The analysis is performed by loading 5 cc of Saltzman reagent into a 
50 cc polypropylene syringe and withdrawing 45 cc of gas sample. The 
syringe is shaken vigorously for approximately three minutes after which 
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Figure 8.1 - Calibration Apparatus 









Percent Transmittance 



Figure 8.3 - Saltzman Calibration Curve 
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time a magenta color completely developed. The contents of the syringe 
are emptied into a cuvette and read on the Bausch and Lomb Spectronic 
20 colorimeter at 550 my against a blank consisting of 5 cc of undeveloped 
Saltzman reagent. Using this technique Meador and Bethea (ref. 19j) report a 
relative error of 0.43% at 41.3 ppm, 1.03% at 26.3 ppm, 2.83% at 9.8 ppm, 
1.77% at 5.7 ppm, and 2.29% at 2.6 ppm. 

Calibration of the Saltzman reagent for N0 2 analysis was done in the 
same manner as for the ammonia analysis. Known concentrations of N0 2 in 
air were prepared by passing air from a compressed air cylinder over a 
Teflon permeation tube filled with liquid N0 2 . The apparatus diagrammed 
in Figure 8.1 was used for this purpose, as it was in the ammonia cali- 
bration work. A typical N0 2 calibration curve is shown in Figure 8.3. 

The Saltzman reagent was prepared by mixing the following chemicals 
with enough distilled and deionized water to make a one liter solution: 

0.050 g N (1-napthyl) ethylene-diamine di hydrochloride 
0.050 g 2-napthol 3,6 disulfonic acid disodium salt 
1.500 g sulfanilamide 
15.0 g tartaric acid. 


8.3 Nitrous Oxide Calibration and Analysis 

Nitrous oxide was analyzed chromatographically using a Varian model 
1532-B trace gas analyzer employing helium ionization detector. The 
following is a listing of the pertinent chromatograph operating condi- 
tion*: which were used: 


packing material Porapak Q 

column length 20 ft 

column temperature 35°C 

detector temperature 100°C 

detector voltage 360 V 

helium flow rate 1 cc/sec 


Because the sensitivity of the helium ionization detectors decayed 
with time, it was necessary to calibrate the chromatograph daily. This 
was easily done with two prepared sample cylinders purchased from the 
Matheson Company containing 33 ppm and 103 ppm N 2 0 in air. Approximate- 
ly four each of these samples were injected daily to provide calibration 
peak heights. Nitrous oxide peak heights from reactor gas samples were 
then referred to these standard peak heights in order to determine N 2 0 
concentration. 

A standard deviation of 0.59 ppm at 33 ppm was calculated from 
typical calibration data. 
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9.0 


NITROGEN OXIDES REDUCTION 


The reduction of nitrogen oxides in the presence of excess oxygen has 
been studied to a considerable extent. Wikstrom and Nobe (ref. 20) studied 
the decomposition of N0 2 over CuO-alumina and a Ce0 2 -alumina catalyst and 
found that even in the presence of oxygen a 45% conversion of 968°F (520°C) 
could be obtained using the CuO-alumina catalyst. Bachman and Taylor (ref. 

21) and Sakaida, et al_. (ref. 22) investigated the high temperature (above 
1000°C) decomposition of NO on platinum, nickel oxide, and platinum oxide 
catalysts. Both invesgi gators found the reaction to be retarded by the pre- 
sence of oxygen. Fraser and Daniels (ref. 23) investigated the decomposition 
of NO over various oxide catalyst at temperatures from 740°C to 1040°C. They 
found the reaction was zero order on all catalysts studied and was very slow 
at temperatures below 800°C. Tamaka and Ozaki (ref. 24) along with Winter 
(ref. 25) have investigated the decomposition of pure N 2 0 over oxide catalysts 
in the temperature range of 400°C to 650°C. Both investigations found the re- 
action to be first order with respect to N 2 0 and to be retarded by the pre- 
sence of oxygen. Amphlett (ref. 26) studied the decomposition of N 2 0 in an 
excess of oxygen and obtained a conversion of 81% at 500°C using a supported 
copper oxide catalyst. Sterbis (ref. 27) studied the decomposition of pure 
N 2 0 on nickel oxide over the temperature range of 300°C to 400°C and found 
the rate was severely reduced by the presence of oxygen. 

All of the work in the literature indicates that for the catalytic mat- 
erials known, the reduction of the oxides of nitrogen (N 2 0, NO, and N0 2 ) in 
the presence of oxygen is extremely slow at temperatures below about 500°C- 
600°C (932°F-1112°F) . At these high temperatures the energy requirements for 
a spacecabin application would be very high. 

The set or reactions by which the removal of the nitrogen oxides might 
proceed is as follows: 


2N0 2 J 2N0 + 0 2 
2N0 J N 2 0 + 3/20 2 
2N 2 0 t 2N 2 + 0 2 


From a theoretical viewpoint, the oxides could be reduced to the thermody- 
namic equilibrium concentration which, since all of the reactions are en- 
dothermic, decreases with decreasing temperature. Using calculated equili- 
brium constants and an atmosphere with an oxygen partial pressure of 0.2 atm 
and a nitrogen partial pressure of 0.48 atm, equilibrium concentrations of 
the oxides were calculated. In order to reduce the N0 2 concentration to 0.2 
ppm, the catalyst would have to be active at a temperature of 370°C (700°F) 
or less. The equilibrium concentration of NO is less than that of N0 2 for 
temperatures below 430°C (800°F). The equilibrium concentration of N 2 0 is 
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extremely low, being less than 10“ 4 5 ppm at 400°C. From the thermodynamic 
calculations it becomes clear that the catalytic removal of nitrogen di- 
oxide down to fractional ppm levels without a reducing reactant puts rather 
severe activity requirements on the catalyst since it must promote the re- 
duction reaction at a relatively low temperature. Unfortunately no catalysts 
are known or were found which are active at temperatures below about 500°C. 

In order to clarify the behavior of nitrous oxide in the presence of the 
ruthenium catalyst at the temperatures expected in a catalytic oxidizer, N 2 0 
was substituted for ammonia at approximately 100 ppm in the reactor inlet. 

At temperatures below 650°F no reaction was detected. Taking into account 
the accuracy of the chromatographic analysis, the N 2 0 disappearance rate could 
not have been higher than 2.6 X 10“ 6 gmole N 2 0/hr/g cat. This rate is not 
high enough to account for either the N0 2 or N 2 which was produced in the am- 
monia oxidation, indicating that N 2 0 is not a precursor for N0 2 or N 2 . This 
finding is in agreement with Schriber's (ref. 5) conclusion that N 2 0 is not re- 
duced to nitrogen at the temperatures investigated. Thus the N 2 0 which is 
formed at these temperatures can be expected to be inert with respect to the 
ruthenium catalyst. 

From the investigations made it appears that the best policy to follow 
for nitrogen oxides control is: 

1. Carefully avoid the formation of these oxides from amnonia 
oxidation by passing all ammonia containing streams through 
an ammonia oxidizer filled with supported ruthenium catalyst. 

This oxidizer should be operated at the minimum acceptable 
temperature which will probably be about 250°F but definitely 
shc-Id •'''•*• hp oDerated above 350°F. 

2. Nitrogen dioxide can be adsorbed on activated carbon with 
some success. 

3. Nitric oxide is rather rapidly oxidized to nitrogen dioxide 
at room temperature and thus should not itself be a problem. 

4. If nitrous oxide formation cannot be avoided, then the only 
known method for its removal in the presence of oxygen is 

reduction at a temperature of about 500°C using either sup- 

ported copper oxide or nickel oxide. 
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10.0 OXIDATION OF AMMONIA IN A STEAM BACKGROUND 


A preliminary investigation was made into the feasibility of carrying 
out the oxidation of small concentrations of ammonia in a stream consist- 
ing mostly of steam. The reactor assembly shown in Figure 10.1 was used 
in this investigation. An ammonia water solution was prepared which con- 
tained 12 to 20 ppm by weight of ammonia. This stream was passed through 
a vaporizer coil and then mixed with air or oxygen, both were used in dif- 
ferent phases of the investigation. The steam-air-ammonia stream then 
passed through a second preheater coil and through the fixed bed reactor. 

After passing through the reactor, the steam was condensed and cooled to 
ambient temperature. During the condensation any ammonia present is ab- 
sorbed back into the water stream. Several analysis of the gas stream 
from the condenser failed to show any ammonia in the gas stream. (De- 
tection limit about 1 ppm of gas by volume). The reactor was always 
operated at atmospheric pressure. 

The catalysts shown in Table 10.1 were investigated for their 
ammonia oxidation efficiency at temperatures up to 200 °C. By far the 
best catalyst of those investigated was the 0.5% platinum on alumina. 

At the temperatures used (140°C to 200°C) the platinum exhibited approxi- 
mately twice the activity of the next best catalyst which was a 0.5% 
palladium on alumina. Both of these catalysts showed a tendency to lose 
activity with use. It is not clear what the mechanism for the catalyst 
deactivation was. The water used in feed preparation was distilled and 
deionized and thus should contain no halogens, heavy metals, or sulfur 
compounds. There is a strong possibility that the water caused a phase 
change in the alumina support thus altering the surface area. Drying the 
catalyst at 200°C did not appear to restore any of the lost activity. 

The oxygen partial pressure did not appear to influence the reaction 
rate in the range from 0.5 atm down to 0.05 atm. 

An extended run was made using a 25 gm charge of 0.5% platinum on 
alumina (Engelhard Industries). The feed to the reactor was a 14 ppm 
solution of ammonia in distilled water and a oxygen flow of 500 standard 
cc/min. The temperature was 146°C (295°F). Under this condition with a 
ammonia-water flow rate of 5.8 ml /min. approximately 97% of the ammonia 
was oxidized to nitrogen and water with no significant formation of nitrate. 

After having processed 13 liters (520 ml solution/gm cat.) of ammonia- 
water solution, the catalyst was oxidizing approximately 94% of the ammonia. 
Therefore the catalyst was losing its activity at a very low rate. The 
effect of impurities in the feed was not determined. The activity at higher 
temperatures was not determined since when the temperature was increased 
above 150°C the ammonia was essentially 100% oxidized. 

Although the investigation was preliminary it does indicate that ammonia 
removal from a steam background is feasible. Catalysts supported on non-alumina 
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TABLE 10.1 

LIST OF CATALYSTS TESTED FOR ACTIVITY IN THE OXIDATION 
OF AMMONIA IN A HIGH WATER BACKGROUND 



CATALYST 

MANUFACTURER 

SURFACE AREA 
m 2 /g 

1. 

0.5% ruthenium on alumina 

Engelhard Industries 

74 

2. 

0.5% platinum on alumina 

Engelhard Industries 

93 

3. 

0.5% palladium on alumina 

Engelhard Industries 

unknown 

4. 

4% nickel oxide, 4% iron 
oxide, 4% cobalt oxide 
supported on alumina 

Harshaw Chemical 

78 

5. 

10% vanadium pentaoxide on 
alumina 

Harshaw Chemical 

115 

6. 

4% manganese dioxide on 
silica 

Harshaw Chemical 

258 
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materials should be invesitgated since it is probable that the deactivation 
experienced may have been due to the support. The deactivation was not so 
severe as to make the process unusable however. Higher temperatures could 
be investigated to determine if undesirable products are produced Also 
the effect of impurities in the feed water should be determined. 
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11.0 CONCLUSIONS 


1. Ammonia in small concentrations (less than 300 ppm) can be oxidized to 
nitrogen and water with essentially no nitrogen oxide formation by 
using a supported ruthenium catalyst and carefully limiting the reactor 
temperature to 250°F or less. At reaction temperatures from 300°F to 
400°F nitrous oxide is a product and at temperatures above 400°F nitro- 
gen dioxide is also produced. 

2. Other good catalysts for the reaction are platinum on alumina and 
Hopcalite. Both of these materials, however, tend to produce more 
oxides of nitrogen as products of the ammonia oxidation than does 
ruthenium. 

3. Quantitative relationships for the ammonia oxidation rate as a function 
of temperature and ammonia concentration are presented. 

4. Ruthenium catalyst does not tend to be deactivated with use as one 
charge operated for over 1400 hours without a loss of activity. The 
reaction rate using ruthenium is not greatly affected by changes in 
humidity level and oxygen concentration. Exposure to hydrogen sulfide 
caused a loss of activity of about 50%. 

5. No catalyst is known which will promote the reduction of nitrogen 
oxides in the presence of oxygen at temperatures below 800°F. Thus 
nitrogen oxide formation must be carefully avoided. 

6. Ammonia in a background of mostly steam can be oxidized to nitrogen 
and water using supported platinum and temperatures of 250°F to 400°F. 
The reaction needs further investigation to determine the effect of 
impurities on reaction rate and catalyst stability. 

7. At 72°F the equilibrium sorption capacity of copper sulfate treated 
sorbeads was significantly greater than that of Hambeads over the range 
of 0-500 parts per million ammonia in the bulk gas phase. 

8. The equilibrium adsorption isotherms can be adequately modelled using 
an expression of the form proposed by Langmuir. These isotherms may 
also be used to estimate the minimum amount of material required to 
adsorb a given load. 

9. The combination of an expression that assumed the surface adsorption 
was controlling with an overall material balance resulted in a mathe- 
matical model which adequately described the observed breakthrough 
curves. A least squares analysis was used to determine recommended 
design values of the rate constant. 

10. One experimental run suggested the possibility of regeneration of 
copper sulfate treated sorbeads. Further work is needed to deter- 
mine the performance of this material under extended cycling. 
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APPENDIX A 


CATALYSTS DATA 

A.l Ruthenium Catalyst Data 


Catalyst Manufacturer Engelhard Industries 

Lot Number 16-988 

Ruthenium Content * 0-5% 

Support Material Alumina 

Catalyst Dimensions 1/8 in. X 1/8 in. 

Surface Area 74 m 2 /9 

Catalyst External Area 10.7 cm 2 /g 

Apparent Bulk Density 60 lb/ft 3 

Total Pore Volume 0-.4 cm 3 /g 

A. 2 Hopcalite Catalyst Data 

Catalyst roanuTacuirei Mine Safety Appliances 

Content Mixed Oxides of Copper and Maganese 

Catalyst Dimension Irregularly Shaped Particles 

Surface Area 156 m 2 /g 

External Area 12 cm 2 /g (Approximate) 
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A. 3 

Catalyst Manufacturer 

Lot Number 

Platinum Content 

Support Material 

Catalyst Dimensions 

Surface Area 

Catalyst External Area .... 

Apparent Bulk Density 

Pore Radius 


Platinum Catalyst Data 

Engelhard Industries 

18,381 

0.5% 

Alumina 

1/8 in. X 1/8 in. 

97 m 2 /9 

10.7 cm 2 /g 

60 lb/ft 3 

25-35 angstroms 
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ADDENDUM 

TO THE FINAL REPORT 


THE CATALYTIC REMOVAL OF AMMONIA AND NITROGEN 
OXIDES FROM SPACECABIN ATMOSPHERES 

By A. J. Gully, R. R. Graham, J. E. Halligan, 
and P. C. Bents en 
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ABSTRACT 


The adsorption of 100 to 500 ppm ammonia from dry air is investi- 
gated experimentally and mathematically. The adsorbent discussed is 
a copper sulfate treated commercial silica gel. The isotherm data is 
well represented by a Langmuir type equation. Using unsteady state 
material balances, a lumped type rate equation, and the Langmuir iso- 
therm, the breakthrough curves are well fitted by adjusting the value 
of the one parameter in the model. A replacement reaction involving 
ammonia substitution for the water of hydration in copper sulfate is 
indicated. 
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Introduction 


The removal of ammonia at concentrations of 500 ppm or lower from 
air is of interest for maintaining life supporting systems such as 
manned space vehicles as well as for the treatment of certain industrial 
vent streams. Several solid sorbents have been investigated for the 
removal of ammonia. These range from silica gel (1), to acid impregnated 
charcoal (2), to chemically treated glass (3). In our laboratory we in- 
vestigated the removal of ammonia from essentially dry air at concentra- 
tions of 100 ppm to 520 ppm using a copper sulfate impregnated silica 
gel as the sorbent. Our objective in the investigation was to sufficient- 
ly correlate the experimental data for the system such that usable pre- 
dictions of adsorber performance could be made. Since adsorbers are 
inherently operated as unsteady state devices, the describing equations 
for the system are partial differential equations in time and 
position. 

Describing Equations 

For a fixed bed adsorber, operating essentially isothermal ly, 
assuming plug flow through the adsorber, negligible pressure drop, con- 
stant velocity, and no radial gradients a material balance for the 
adsorbed component of interest yields equation (1). 

_v If = R + e B It (i) 

3 

where C = gas phase adsorbate concentration mole/cm 
Z = bed position, cm 
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,T = time, min 

v = superficial bed velocity cm/mi n 

R = rate of adsorption, mole/cnr-ymin 

e B = bed external void fraction. 

A material balance for the adsorbate loading on the solid yields 
the defining equation for the rate of adsorption. 

P B 3lf = R ^ 

3 

where p g = bulk density of the packed bed, g/cm 
W = adsorbent loading, mole/g. 

In a detailed, mechanistic analysis the rate of adsorption on the 
surface of the solid is no doubt an involved function of the fraction of 
the surface covered. At a given temperature and gas phase composition 
there is a fixed amount of adsorbate that can be held on the solid sur- 
face. It is reasonable, then, to say that one measure of the driving 
force for surface adsorption is the local displacement from equilibrium. 
The rate of adsorption can often vary with position within an individual 
adsorbent particle due to the existence of concentration gradients within 
a particle. However in this study the intent was to use as manageable 
a describing model as possible. The first attempt thus was to use a 
lumped type equation to describe the rate of adsorption. The rate 
equation used was 


R = k (W E -W) 


( 3 ) 


3 

where k = adsorption rate constant g/cm /min 

= equilibrium loading, based on local gas phase adsorbate 
concentration, mole/g . 
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Several models have been proposed to relate the equilibrium solid 
loading to the gas phase adsorbate concentration at constant temperature. 
Freundlich and Langmuir isotherm forms were considered for use in this 
work. The Langmuir form was found to yield a superior fit of the data 
and was used in the further analysis of the adsorber. Thus the equili- 
brium equation used was 


W 


E 


aC 

1 + bC 


(4) 


where a and b are temperature dependent constants. 

For the work reported here the inlet concentration to the adsorber 
bed was stepped at time equal zero from a zero ammonia in air concentra- 
tion to some preset constant concentration value. The adsorbent bed 
was in all cases but one an unexposed sample of copper sulfate treated 
silica gel. Thus the system initial and boundary conditions were 


t r\ +. \ 

= r 

in 

(5) 

C(Z,0) 

= 0 

(6) 

W(z,0) 

= 0 . 

(7) 


The approach in this work was then to simulate the transient behavior 
of the adsorber using these seven equations which contain seven system 
constants. Values for all seven constants must be chosen a priori to a 
numeric solution of the equation set. A preferred approach might be to 
reformulate the problem such that the variables are dimensionless and such 
that a minimum number of parameters are involved in the problem solution. 
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The following changes of variables accomplish these objectives: 


t = tk/p B (8) 


Z = ZW F k/v/C. n (9) 


C = c/c in (10 


W = w/w F (11 


Wp = W E /Wp . (12 


After making the above substitutions the following two parameter equation 
set results: 


8C _ 9W , e B C in, dC 

'' ^ + \ n w ) * 

az at p b f at 


(13 


/N A 



at l + bc in c 


c(o,t) = 1 
c(z,o) = o 

W(Z,0) = 0 . 


(14 

(15 

(16 


(17 


5 


For the physical system under consideration the last term in the gas 
phase material balance was found to be insignificant since 

6 * ^ 4 ^ * '°' 6 • (18) 

If the gas phase accumulation term is dropped the remaining equation set 
contains only one parameter, bC in - For the special case where bC in equal 
zero (a linear isotherm), the linear equation set can be Laplace trans- 

A A 

formed with respect to t and integrated with respect to Z to give the 
following equation in the Laplace domain: 

LC = exp [-Z (^)Vs . (19) 

By using the series expansion for the exponential 

2 3 

e X = l+ X+ |j + ||- + etc. (20) 


the following time domain solution for the linear case can be found 

t 


C = 1 - 


oo oo 



(-1)^1 Z n+ m t m e -t 
(m!) 2 (n-1 ) ! (n+m) 


( 21 ) 


This analytic solution was utilized to help check the stability of the 
numerical solution program. 
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Adsorber Loading at Bed Inlet 

The inlet of the adsorber is exposed to a constant gas phase ammonia 
concentration. Thus the loading as a function of time can be obtained 
analytically. At the inlet equal one (note that this result is inde- 
pendent of isotherm form) thus 


1 - i 

dt 

(22) 

A 

W(at Z = 0) = l - e _t . 

(23) 


For an adsorber operated as described here, the bed loading is most rapid 
at the inlet, therefore any external mass transfer resistances would be 
most significant at the beginning of the cycle. 

Numerical Solution Technique 

Because of the non-linear nature of the problem, numerical solution 
techniques were used to solve the dimensionless equation set for various 
values of the parameters 3 and C^ n . The approach used was to apply a 
Euler type finite difference equation to start the solution and then switch 
to second order finite difference formulae for the remainder of the calcu- 
lations. For the relatively short beds and low concentrations used in 
this study, the gas phase accumulation term was found to be totally 
insignificant and was dropped. 

Comparison with Experimental Data 
The experimental data collected was in the form of concentration 
history (or breakthrough) curves. These curves are a plot of the bed exit 


A8 



ammonia concentration as a function of run time. Calculated bed capa- 
cities were determined from the integrated average exit ammonia con- 
centration. The isotherm data thus generated are shown in Figure 1 
for the two temperatures investigated. The data are plotted a C vs 
C/W which for the Langmuir form should result in a linear relationship 
of intercept -1/b. In order to compare the experimental data to the 
numerical solution, the unknown system constant k had to be determined. 

A /V 

Since for the particular dimensionless variables chosen both t and Z 

A 

contain k, the numerical solution was utilized as a plot of C versus 

A A A 

t/Z with lines of constant Z. Typical numerical solution results are 
shown in Figure 2. As may be seen from this figure, the solution is 
not greatly sensitive to changes in the parameter bC in in the range from 
bC in equal about 2 to bC.- n equal about 4. 

The central problem was to determine the value of k which would 
best represent the entire data set at each temperature. A crude esti- 
mate of k can be made using the initial breakthrougn ammonia uum.cn- 
tration. For the early minutes of the run, the bed loading is essential- 
ly zero for all lengths. Thus the material balance can be integrated 
assuming W = 0. 


3C 

“ x- 

3Z 


„ 1 + bC. 

- C [ 


1 + 


bC. C 
in 


(24) 


Z = 


Z Wp k 


1 


v C 


in 


T+ bC 


— CbC in (l - C) - In C] 


in 


(25) 
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' FIG. l - LANGMUIR ISOTHERM FIT 
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FIG. 2 - effect of isotherm parameter 



The algorithm used to obtain the best value of k is outlined below: 

(1) The value of k was crudely estimated from the initial break- 
through ammonia concentration by the use of equation (25). 

(2) The equation set was numerically solved for each value of the 
parameter bC^ and the results retained for values of k from 
one-half the smallest estimate to twice the largest estimate. 

(3) A trial and error least squares fitting program was used in 

A 

which a value of k was selected, the value of Z calculated for 
each run, and the summation of deviation squared determined. 

(4) Step (3) was repeated until the minimum sum of squares was 
determi ned . 

Using this least squares procedure the values of k determined were 

k = 0.0053 — j 2 — at 72°F 
cnr mi n 

k = 0.0214 g/cm 3 /min at 100°F . 

The major experimental conditions for the adsorption runs considered are 
shown in Table 1. In addition the equilibrium loading and value for the 
parameter bC in are shown. Figures 3 through 8 show the calculated break- 
through curves with data points superimposed. Since the solution is only 
weakly dependent on the value of the parameter bC in , the calculated break- 
through curves for several of the runs are essentially the same, and these 
are shown on a common plot. The agreement between predicted and experimental 
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TABLE 1 - EXPERIMENTAL CONDITIONS 


Number 

V 

cm/mi n 

Z 

cm 

W F 

g/g 

C in 

ppm 

P=BC i 

21 

1340 

3.15 

0.029 

290 

2.13 

22 

1365 

1.92 

0.021 

150 

1.10 

24 

1395 

3.91 

0.027 

325 

2.39 

25 

1509 

1.98 

0.030 

520 

3.82 

26 

1314 

2.19 

0.016 

90 

0.66 

27 

1320 

1.91 

0.022 

215 

1.58 

28 

1269 

2.04 

0.009 

235 

.74 

30 

1388 

2.46 

0.012 

440 

1.39 

32 

1223 

2.37 

0.008 

210 

.66 
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FIG. 3 - COMPARISON OF EXPERIMENTAL DATA AND 
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CALCULATED BREAKTHROUGH 
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curves is reasonably good for all runs except number 25. The cause for 
the poor fit for run 25 is not clear. This run had the highest inlet 
ammonia concentration. The point for run number 25 is well fitted by 
the isotherm equation which tends to rule out experimental error in the 
flow measurement, solid charge weight, and average exit concentration. 

For reference the breakthrough curve is given for a rate of adsorption 
twice the least squares value. For the 72°F data as a whole, a five 
per cent variation ink had an insignificant effect on the least squares 
calculation. 

A color interface (aqua to dark blue) traversed the bed early in 
each run. Bed exit ammonia concentrations were less than 20 per cent of 
the inlet concentration during the traverse. There was, however, no 
detectable discontinuity in the breakthrough curves as the color inter- 
face exited the bed. The color change was thought to be due to the 
replacement reaction (4) 

CuS0 4 • 5H 2 0 + 4NH 3 t CuS0 4 • 4NH 3 + 5H 2 0 . 

The characteristic color of CuS0 4 • 4NH 3 is dark blue. Another indication 
that some type replacement occurred during a run was that the bed typical- 
ly gained only 86% to 93% of the weight of airmonia adsorbed. 

The latent heat of adsorption was determined using the two isotherms 
to be approximately 22 kcal/mole. Since the latent heat of condensation 
of ammonia is only about 6.5 kcal/mole, an activated adsorption is again 
indicated. 
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Regeneration Study 

Only a very limited investigation was made on the regeneration of 
saturated sorbent. Adsorbent from run 24 was regenerated by heating 
in a cross flow oven for several hours at a temperature of 385°F. The 
sorbent was then resaturated at identical conditions to run 24. The 
adsorbent capacity after one regeneration was 80% of the fresh adsor- 
bent's capacity. The characteristics of the breakthrough curve for 
the regenerated sorbent were similar to the fresh sorbent except with 
a slightly lower rate. 

Conclusions 

Copper sulfate impregnated silica -gel compares favorably with other 
ammonia adsorbents. Acid impregnated charcoal is reported to have an 
ammonia capacity of 2.5% (mass) at a bulk ammonia concentration of 
moo Dorn. The treated silica gel used in this study had a capacity of 
3% (mass) at 520 ppm ammonia. 

The isothermal breakthrough curves for the relatively short. beds 
used can be well represented using a single parameter equation set. 

The model used assumes surface adsorption is the rate controlling step 
and that mass transfer resistances are minor. By putting the equation 
set in dimensionless terms only a single parameter remains, thus making 
the numerical solution more general. 

The copper sulfate impregnated silica gel can be regenerated to 80% 
of its original capacity by heating at 385°F. 
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